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1.  INTRODl'CTION 

Wave  drag  reduction  is  one  of  the  imi)ortant  goals  of  designing  a high  perfor- 
mance aircraft  capable  of  operating  in  the  supersonic  regime.  Since  there  are  many 
design  requirements  which  conflict  with  one  another  as  far  as  wave  drag  reduction  is 
concerned,  an  optimization  procedure  is  needed  to  determine  the  minimum  wave 
drag  configuration  subject  to  the  constraints  imposed  by  these  requirements.  Until 
recently,  the  procedure  relied  heavily  on  experience  gained  through  actensive  wind 
tunnel  testing  of  various  geometries.  Such  a design  procedure,  which  usually  could 
not  be  carried  out  systematically,  was  expensive  and  time  consuming.  However, 
the  advancec  of  recent  years  in  numerical  methods  and  computer  technology  have 
made  feasible  systematic  optimization  procedures  using  ecact  numerical  methods. 
Consequently,  a wave  drag  reduction  proc«iure  using  the  method  of  characteristics 
has  been  developed  which  is  presented  in  this  report. 

The  present  wave  drag  reduction  procedure  makes  use  of  t\vo  basic  methods: 
the  Latin  Square  sampling  technique  and  the  Three-Dimensional  Method  of  Character- 
istics. The  former  is  used  to  select  sample  configtirations  so  efficiently  that  a 
small  number  of  samples  "san  well  represent  the  entire  family  of  configurations. 

The  latter  is  used  to  calculate  accurately  the  wave  drags  of  the  sampled  configura- 
tions, Briefly  stated,  the  present  approach  consists  of  calculating  the  wave  drag 
of  a baseline  configuration  and  some  variations  specified  by  the  Lritin  Square  sampl- 
ing technique,  determining  a functional  dependraice  of  the  wave  drag  on  these  varia- 
tion parameters,  and  minimizing  this  wave  drag  function  to  obtain  the  configuration 
with  minimum  wave  drag.  The  procedure  is  general  with  respect  to  the  number  of 
geometric  parameters  (or  variables);  the  higher  the  number,  the  larger  the  required 
Latin-Square  size.  The  computer  programs  developed  under  this  study  cover  the 
mr)st  often  used  3x3  and  5x5  Latin  Squares. 
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The  complete  wave  drag  reduction  prograin  has  been  carried  out  in  two  phases. 

In  the  fii'st  plnise,  a i^rocedure  was  developed  for  minimizing  the  wave  drag  of  a 
forward  fuselage  and  canopy  configuration  as  reported  in  Reference  1.  In  the 
stx’ond  phase,  the  procedure  has  been  aKp;uided  to  account  for  the  influence  (jf  the 
wing  :md  wing-body  blending  on  the  overall  wave  drag.  In  this  final  report,  the 
rcseai’ch  performed  under  the  wave  drag  reduction  contract  is  presented.  The 
phase  I work,  which  was  reported  in  Reference  1,  i.s  included  as  Appendix  A for 
rcitdy  reference,  whereas  the  main  text  of  the  report  presents  the  phase  II  work 
and  touches  upon  some  of  the  phase  I work.  A 3 x :i  Latin  Square  was  used  in  phase 
II  ;ind  a 5 X 5 Latin  Square  in  phase  I . The  surface  fitting  method  using  Latin 
Squares  as  presented  in  Reference  2 was  improved  during  the  phase  I study:  the 
improvement  which  is  presented  in  Appendix  A,  is  essential  for  the  success  of  the 
optimization  procedure  (here  applied  to  wave  drag  reduction).  The  procedure  using 
the  improved  Latin  Square  surface  fitting  method  has  been  provai  in  both  phases 
through  application  to  the  F-4  configuration.  For  further  validation  of  the  procedure, 
it  was  applied  to  the  von  Karman  ogive.  For  given  configuration  length  and  base, 
the  present  optimization  procedure  correctly  predicted  the  von  Karman  ogive  as 
the  minimum  wave  drag  body. 

In  this  report,  the  basic  approach  is  given  first,  which  consi.sts  of  the  formu- 
lation of  the  problem  and  a brief  account  of  the  two  basic  methods.  It  is  followed  by 
a discussion  of  the  method  of  describing  the  body  and  the  selection  of  geometric 
variables  and  their  ranges  for  defining  a family  of  configurations.  Then  the  flow 
field  calculation  and  the  wave  drag  equation  are  presented.  Sample  residts  of  the 
calculated  flow  fields  are  given,  and  calculated  wave  drag  coefficients  are  trdjulated. 

These  coefficients  were  used  to  derive  the  wave  drag  equation  which  expresses  the 
wave  drag  as  a function  of  the  geometric  variables.  Once  the  wave  drag  equation  is 
obtained,  the  dependence  of  the  wave  drag  on  the  geometric  variables  is  established. 
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A set  of  fi(^\u-es  are  {pivcn  to  illustrate  some  of  the  characteristics  of  the  wave  dra;; 
■equation.  This  is  followed  by  a presentation  of  the  optimization  procedure  and  the 
prediction  of  the  minimum  wave  drag  body.  Finally,  conclusions  are  drawn  and  some 
recommendations  given.  The  Latin  Sejuare  technique  including  the  method  of  con- 
struction is  presented  in  Appendbe  B.  A discussion  of  the  general  body  description 
method  is  given  in  Appendix  C.  The  Numerical  Search  Procedure  for  the  minimum 
wave  drag  configuration  is  presented  in  Appendix  D.  The  validation  of  the  optimiza- 
tion procedure  using  the  von  Karman  ogive  is  presented  in  Appendix  E. 
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-•  1 oltMl'l-VrioX  AM)  METHODS 

The  wave  di'U'j;  >'l  a coiit'is^uralion  is  a i'unctinn  nt  ;•  number  of  factors.  For 
Hivon  flight  conditions  tho  wave  dra^  di*j»ends  on  how  the  eonfi^ration  is  shaped. 
Since  the  sliapo  of  an  airei  ait,  at  !e;ist  in  the  initial  desi^^n  phase,  is  mainly 
iletermined  by  considerations  other  than  the  wave  ora;;,  it  is  practical  to  consider 
the  im)l)lcm  of  reducing;  the  wave  di-ap:  of  a <;iven  baseline  confiv-i  ation  by 
obtainin',;  a v-ariation  confi^oiration  thai  satisfies  all  the  dcsi^o  const  naints  yet  h:is 
the  least  wave  dra;;.  Such  a baseline  confi;;uration  could  be  a new  confiiniration 
at  a cei'tain  stage  of  devebipment  <>r  it  could  be  an  evisting  airplane  that  is  to  be 
nKxiifiod  or  improved. 

In  this  section,  the  formulation  of  the  wave  drag  reduction  problem  is  outlined, 
and  the  two  basic  methods  to  be  used  in  this  study  are  introduced. 

a.  Formulation 

The  baseline  c<»ifiguration  can  be  described  by  a set  of  geometric  variables. 

A family  of  configurations  including  the  baseline  can  be  generated  by  assigning 
different  values  to  some  or  ail  of  the  geometric  variables.  If  the  wave  drag  can  be 
expressed  as  a function  of  these  variables,  a particular  set  of  values  of  these 
variables  that  gives  the  least  wave  drag  can  be  found  by  minimizing  the  function. 

This  set  of  values  then  ])roduces  the  minimum  wave  drag  configuration. 

The  key  to  this  problem  is  how  to  obtain  such  a functional  cK))ression  for  the 
wave  drag.  For  the  present  study,  four  geometric  variables  are  considered  in 
defining  the  family  of  configurations.  If  each  varialrle  assumes  three  values,  the 
evaluation  of  the  partial  derivatives  with,  respect  to  these  four  variables  for  a Taylor- 
series  type  expression  would  require  81  wave  drag  calculations.  It  will  be  shown 
that  through  Latin  Square  sampling  the  presatt  procedure  proves  useful  using  only 
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10  wave  drag  calculations.  In  the  present  aj)proach,  the  wave  drag  coefficient 
is  first  assumed  to  be  of  the  form 

4 

^DW  ^ ^o  ' ( Vi  ' ) ‘ (1) 

1 1 ' 

where  x.(i  1,4)  are  the  geometric  variables  and  a.  and  a.,  are  to  be  determined. 

1 ’ " i ij 

o 

The  Latin  Scjuare  sampling  technique”  is  used  to  sample  0 sets  of  v'alues  of  x. 
out  of  the  total  population  of  81  sets.  The  wave  drag  coefficients  for  the  F-4  type 
baseline  and  the  configurations  defined  by  each  of  the  9 sets  arc  then  calculated 
by  the  Three-Dimensional  Method  of  Characteristics. ' ’ When  the  10  calculated 
wave  drag  coefficients  and  the  corresponding  /alues  of  the  geometric 

variables  a.  are  substituted  into  Equation  1,  n linear  equations  for  the  10  un- 
knowns a^  and  are  obtained.  These  equations  arc  then  solved  for  the  a.  and 
Ujj,  which  are  substituted  back  into  Equation  1 to  produce  the  functional  axpres- 
sion  for  the  wave  drag  in  terms  of  the  geometric  variables.  By  minimizing 
in  Equation  1 subject  to  a given  set  of  cr  »straints,  e.  g. , a given  volume 
of  the  aircraft,  the  minimum  wave  drag  configuration  corrcs]>onding  to  the  given 
set  of  constraints  is  'etermined.  In  the  present  study,  a numerical  search  pro- 
cedure is  used  to  tind  the  minim*vm  wave  drag  configuration  (Appendix  D), 

b.  I.atin  Scuare  Sampling  Technique 

The  Latin  Square  method,  which  has  mostly  been  used  in  agriculture  and 
biological  research,  is  a very  efficient  sampling  technique  and  is  much  better  than 
random  sampling.*  For  this  study  a particular  type  of  Latin  Square  (the  orthogonal 
squares)  suitable  for  a variety  of  technical  problems^  is  adopted.  With  this  tjq)e 
of  Latin  Square  arrangement,  a .1  x 3 square  is  the  correct  size  for  four  geometric 
variables  x^,  each  taking  three  values.  It  is  convenient  to  introduce  the  reduced 
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variables  z.,  which  are  related  to  the  geometric  variable  x.  through 

X.  X.  z.  - X.  \ 

1 1 . Ill  1 . I 

max  min  ^ ' max  min ' 

■'^i  2 " 2 ~ 

where  subscripts  max  and  min  denote  the  maximum  and  minimum  values, 
respectively.  Corresponding  to  the  three  values  of  the  geometric  variables  x., 
the  reduced  variables  z.  always  aseume  the  three  levels  of  0,  11.  The  3x3 
Latin  Square  arrangement  in  terms  of  the  levels  of  the  reduced  variables  and 
the  cell  number  is  shown  in  Table  1.  * It  is  seen  that  in  this  way,  the  Latin 
Square  arrangement  remains  the  same  whatever  the  values  of  the  geometric 
variabtes. 

c.  Three-Dimensional  Method  of  Characteristics 

In  this  study  the  Three-Dimensional  Method  of  Characteristics^’^  is  used  to 
calc-ilate  the  wave  drag  coefficients  for  the  fuselage-wing  configurations  sampled 
by  the  Latin  Square  technique.  This  method  has  been  previously  applied  to  cal- 
culate the  flow  fields  over  a wide  var^'ty  of  configurations  including  spherically- 
capped  three-dimensional  bodies  and  wings^,  aircraft  fuselages  and  wings  at 

g 

general  angles  of  attack,  and  slab  delta  wings  and  space  shuttle  wing-body  con- 
7-9 

figurations.  Whenever  experimental  data  were  available  for  comparison,  good 
agreement  between  theory  and  experiments  was  observed.  The  capability  of 
treating  the  canopy,  however,  was  developed  during  phase  I of  this  study. 

Tn  a related  research  program,  thcThrce-Dimcnsional  Method  of  Characteristics 
was  further  extended  and  improved  to  trctU  realistic  aircraft  wing-body  configura- 
tions including  wing-body  blending.  With  some  modification  and  adaptation,  this 
improved  method  was  applied  to  calculate  the  flow  fields  and  wave  drags  of  the 
variation  configurations. 


*A  discussion  of  the  Latin  Square  construction  is  given  in  Appendix  B- 
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3.  DESCRIPTION  AND  VARIATION  OF  FUSELAGE- WING  CONFIGURATION 


Three-dimensional  bcxiy  description  requires  a great  deal  of  effort  which  at 
times  becomes  extremely  tedious.  Basically,  two  types  of  description  can  be 
made,  analytical  and  numerical.  The  latter  can  describe  complicated  geometry 
accurately,  but  is  cumbersome  for  preparing  input  data.  The  former  is  much  sim- 
pler to  input  and  can  define  a family  of  configurations  based  on  a few  geometric 
variables.  Hence,  analytical  description  is  chosen  for  the  present  study. 

a.  Body  Description 

Every  configuration  has  a number  of  generating  lines,  such  as  the  upper 
profile,  the  lower  profile,  the  maximum  breadth  line,  or  the  wing  leading  edge. 

In  the  present  body  description  procedure,  each  generating  line  is  divided  into  a 
number  of  segments  to  permit  each  segm«tt  to  be  described  by  a conic-section 
curve.  At  each  cross  section  of  the  configuration,  simple  analytic  curves,  e. g. , 
the  ellipse  or  cubic,  connect  any  two  adjacent  generating  lines  to  form  the  contour 
of  the  cross  section.  The  configuration  is  thus  described  analytically  by  simple 
low-order  curves.  For  a smooth  body  a unique  normal  to  the  surface  exists  every- 
where, and  this  condition  usually  requires  slope  continuity  at  the  junctures  between 
two  contour  curves  or  two  segements  of  a generating  line. 

The  fuselage  is  located  in  a right-handed  coordinate  system  where  the  Y-axis 
is  aligned  with  the  fuselage  axis;  the  X-axis  is  spanwise  and  the  Z-axis  is  up.  All 
generating  lines  are  represented  by  a general  curve  fit  of  conic  sections  in  several 
segments.  The  conic-section  curve  takes  the  form 
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be  divided  into  as  many  se  gments  as  necessary  to  provide  adequate  body  description. 
Each  segment  must  be  continuous  with  the  previous  segment  and  with  very  few 
exceptions  the  slope  must  be  continuous  at  the  junctures  to  satisfy  the  requirement 
of  a unique  normal  to  the  surface. 

A typeial  cross  section  of  the  wing-body  configuration  is  shown  schematically 
in  Figure  1.  The  contour  of  the  cross  section  begins  with  a straight  line  represent- 
ing the  canopy  flat.  The  canopy  contour  from  points  C to  0 is  circular  but  can  be 
elliptic  in  general.  The  upper  fuselage  is  represented  by  an  elliptic  curve  from 
points  U to  M.  A straight  side  flat  from  points  M to  Fj  joins  the  upper  fuselage 
to  the  lower  fuselage,  which  is  also  represented  by  an  elliptic  cuiwe  from  F^  to  L. 

A straight  line  from  point  L to  the  centerline  describes  the  bottom  flat.  For  the 
description  of  the  wing,  straight  lines  E G and  F H represent  the  upper  and  lower 
surfaces  of  the  wing,  respectively.  Partial  ellipses  G I and  H I complete  the  wing 
description  near  the  leading  edge  I.  The  wing-body  blending  is  effected  by  circular 
arcs  34  and  56  with  radii  r^  and  rj^,  respectively.  A further  discussion  of  the 
body  description  method  illustrated  by  the  description  of  the  fuselage  and  canopy 
of  Phase  I Is  presented  in  Appendix  B. 

b.  Fuselage  and  Canopy 

The  equation  of  the  canopy  is  given  by 


z-z^oo 

2 

x-x^oo 

z^oo  - Z^(Y) 

X^(Y)  - X^(Y) 

2 2 

If  it  is  a circular  arc,  then  (2^-  Z ) « (X  X , The  equation  for  the  upper 

^ . C , V >' 

fuselage  is 


Z - Z^OO 

2 

1 

■ X - x^joo 

Zy(Y)  - Zj^(Y) 

r 

x„(Y)  - 

1 = 0 


(5) 


The  intersection  between  the  canopy  and  fuselage  is  faired  by  a cubic  from 
points  1 to  2 (Figure  1).  The  projections  Z^^fY)  and  of  lines  1 and  2 on  the 
Y-Z  plane  are  given  by  Equation  3 where  the  coefficitaits  P,  Q,  R,  S,  T are  input 
quantities.  The  projections  (Y)  and  X^fY)  on  the  X-Y  plane  are  obtained  by 
solving  Equations  4 and  5,  respectively. 


X^(Y) 


X^fY) 


(6) 


(7) 


The  fairing  curve  matches  the  slopes  of  the  ellipses  at  both  end  points  1 and  2.  The 
slopes  are  obtained  by  differentiation  of  Equation  4 and  5 


X’ 


/ ax  \ 

I^C 

1 a z /, 

i Ui-^c‘ 

)■ 


The  cubic  equation  that  satisfies  Equations  6 and  8 can  be  written 


(8) 


(9) 


(10) 


The  coefficients  c and  d are  obtained  by  applying  Equations  7 and  9 


c » 3 (X2  - Xj)  - (X^  + 2Xp  (Zg  - 2j)  (11) 

d = -2(X2  - Xj)  + (X^  ^ Xp  (Z^  - Zj)  (12) 

Equation  10  with  c and  d given  by  Equations  11  and  12  is  then  the  cubic  equation  for 
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the  fairing  curve  from  points  1 to  2.  The  quantities  X^,  XJ^  and  X^  are  given 
by  Equations  6 and  7 while  X^,  X^j.  X^,  X^,  Zy,  Zj^,  Z^,  Z^,  Z^,  and  Z^ 
are  obtained  from  Equation  3 where  corfficients  P,  Q,  R,  S and  T are  input 
quantities. 

The  lower  fuselage  is  descnbed  by  an  equation  similar  to  Equation  5 

2 

- 1 - 0 (13) 


’ z - z 

^f 

2 

X-Xj^(Y) 

Zj  (Y)-Z„(Y) 

la  * £ 

Xe^OO  - X^CO 

All  flats  are  given  by  simple  straight-line  equations, 
c.  Wing  and  Blending 

The  upper  and  lower  surfaces  of  the  wing  are  given  by  straight-line  equations 
for  E G and  F H.  Near  the  leading  edge  I,  the  partial  ellipses  G I and  H I are 
derived  as  follows.  The  equation  may  be  writtKi  in  the  form 


(X  - Xj  > a)2  (Z  - Zj)^ 
a b 


(14) 


where  a and  b are  the  axes  to  be  determined.  Differentiation  of  Equation  14  gives 


X - Xj  + a 


Z - Z 


I d^ 

dX 


(IS) 


at  point  G,  X * Xq,  Z = Zq,  and  slope  of  line  E G . Sub- 

stituting these  values  into  Equation  14  and  15  leads  to 


(Zq  - _ 


(14a) 


\dxL 


* 0 


(15a) 
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Elimination  of  from  Equations  14a  and  14b  yields 


2 

Elimination  of  a from  Equations  14a  and  14b  leads  to 

(IT) 

where  a is  given  by  Equation  16.  Equatitm  14  with  a and  b givoj  by  Equations  16  and  17 
describes  the  partial  ellipse  G I.  By  changing  the  subscript  G to  H,  the  equation  for 
the  partial  ellipse  H I is  obtained. 

The  projections  of  lines  passing  through  E,  F,  G,  H and  I on  the  Z-Y  plane 
and  X-Y  plane  are  input  quantities  for  defining  the  wing.  In  order  for  the  partial  ellipses 
to  exist,  points  G and  H must  be  located  within  a certain  range,  which  depends  on  the 
relative  positions  of  these  five  points.  Wh«i  the  wing  span  is  very  small,  it  is  difficult 
to  input  both  projections  of  lines  through  G and  H such  that  these  points  are  located  within 
acceptable  ranges.  In  such  cases,  the  X coordinates  of  G and  H are  calculated  inside 
the  program  to  satisfy  the  range  requirement. 

The  upper  and  lower  blmiings  between  the  fuselage  and  the  wing  are  described 
by  circular  arcs  with  radii  r^  and  rj^,  respectively.  When  the  blending  radii  r^  and 
are  specified  by  geometric  variables  (see  Section  3d.),  points  3 and  5 of  Figure  1 
can  be  obtained  numerically  through  an  iteration  procedure.  However,  in  order  to 
obtain  analytical  normals  to  the  blending  surface,  an  anal3rtlc  expression  must  be  derived. 
Hence,  the  following  procedure  was  used  and  is  illustrated  by  the  upper  blending.  The 
numerically  obtained  line  passing  through  point  3 is  considered  a generating  line.  The 
projection  (Y)  on  the  Y - Z plane  is  expressed  in  the  form  of  Equation  3,  where 
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the  coefficients  are  input  quantities.  At  a given  fusDiage  station,  X^OO  are  obtained 


by  solving  Equation  5. 


X3OO  = (X^-Xy)  ) 


2 11/2 


U M 


The  slope  of  the  tangent  at  point  3 can  also  be  obtained  from  Equation  5 as 


^ Si  - ""u 


) m,) 


Thus,  the  equation  of  the  tangent  3P  (see  Figure  1)  at  point  3 is  given  by 


Z . MjOt-Xj)  .Zj 


The  equation  of  line  E G can  be  written  as 


Z = M^(X  - Xg)  + Zg 

Z„-  Zp 

where  M . = ~ 

^ “ -^E  I 

Intersection  P of  these  two  lines  is  bbtained  from  Equations  20  and  21 


X » 
P 


"^3  * ^4  ! 


Z « 
P 


M„  - 
3 4 


The  equation  of  the  line  bisecting  the  angle  / 3PG  is  given  Ir* 


Z=*M.  (X-X)+Z  ! 

5 ' p'  P I 


tan”^  Mg  + taii  ^ ^^4  j 


where  M, 


and  M is  ^iven  by  Equation  19  and  M by  Equation  22.  The  intersection  between 
the  bisecting  line  and  the  normal  at  point  3 

X - Xg  * M^(Z  - Z3)  = 0 (25) 

is  the  center  0 of  the  arc;  hence  from  Equations  24  and  25 


^0  = 


1 . M3M3 


® * 1 * M^Mg 


(26) 


where  X and  Z are  given  by  Equation  23,  M«  by  Equation  19  and  M.  by  Equation 

p p a V 

24a.  The  equation  of  the  circular  arc  for  the  upper  blending  is  thus 


(X-x/m2-Z„)2  - <27) 

where  X^  and  are  given  by  Equation  26,  Notice  that  the  radius  is  slightly 
different  from  r^  because  of  slight  errors  introduced  by  the  iteration  procedure 
in  obtaining  point  3 and  by  the  fitting  of  the  generating  line  passing  through  that  point. 

d.  Geometric  Variables  and  Variations 

Four  geometric  variables  were  selected  to  generate  a family  of  fuselage-wing 
configurations  including  the  F-4  baseline.  As  illustrated  in  Figure  2,  these 
variables  are  the  horizontal  displacement  a of  the  maximum  horizontal  breadth 
line,  the  lower  deck  height  h (which  itnreases  as  the  lower  profile  is  raised), 
and  the  blending  radii  r^^  and  at  F.  S.  280  and  F.  S.  360,  respectively.  All 
configurations  generated  by  these  variables  have  the  same  canopy  as  the  F-4  and 
must  satisfy  the  over-the-side  view  line  limitation.  The  correspondence  between 
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the  geometric  variables  and  the  reduced  variables  defined  in  Equation  2,  the 
ranges  of  variation  of  these  geometric  variables,  and  the  baseline  values  of  these 
variables  are  tabulated  in  Tab’e  2. 

The  geometric  variables  specify  the  blending  radii  at  two  fuselage  stations 
only.  In  order  to  fully  describe  the  blending  between  the  fuselage  and  wing,  a pro- 
cedure is  needed  to  provide  the  blending  radius  systematically  at  any  given  fuselage 
station  in  the  blending  region.  The  blending  radius  distribution  as  a function  of  the 
fuselage  station  must  satisfy  the  following  conditions: 

1.  At  fuselage  stations  220  or  430,  the  radius  is  equal  to  the  minimum  radius  of 
the  baseline. 

2.  At  F.S.  280  and  F.  S.  360,  they  are  equal  to  r^  and  respectively. 

3.  Bet\veen  F.  S.  280  and  F.  S.  360,  the  radius  must  not  overshoot:  i.  e. , it  must 
not  exceed  the  greater  of  r^  and 

4.  When  rj^  and  are  equal  to  the  minimum  radius,  the  baseline  must  be 
recovered;  i.  e. , the  radius  must  be  constant  throughout. 

5.  When  either  r^^  or  rg  is  equal  to  the  minimum  radius,  the  blending  radius  must 
not  undershoot  anywhere. 

Figure  3 illustrates  some  of  the  possible  radius  distributions  and  serves  as  a 
reference  for  the  following  discussion  of  the  procedure.  The  radius  distribution  is 
given  by  three  cubic  interpolation  formulas  for  the  three  intervals.  The  slope  of  the 
distribution  curve  is  zero  at  F.  S.  220,  280  and  360  for  r » 1. 5 and  at  F.  S.  280  and 
360  for  r = 61. 5.  At  F.  S.  280  the  slope  at  r » 46. 5 is  assumed  to  be  equal  to  that 
of  the  straight  line  At  any  other  r the  slope  is  obtained  from  a spline  fit  of  the 
slope  versus  r in  such  a way  that  the  spline  curve  passing  through  these  three  points 
with  assumed  zero  curvature  at  both  end  points.  The  slopes  at  F.S.  360  are  obtained 
in  an  anologous  way.  Thus,  for  a given  pair  of  r^  and  r^,  the  radius  distribution  curve 
is  determined  by  three  cubic  interpolation  formulas:  the  first  cubic  passes  through 
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r ^ 1. 5 at  F.  S.  220  with  a zero  slope  and  r - at  F.  S.  280  with  a slope  Riven 
by  the  spline  fit;  the  second  cubic  passes  throuRh  r ^ at  F.  S.  280  with  the  same 
slope  as  the  first  cubic  and  through  r = r2  at  F.  S.  300  with  a slope  given  by  the 
second  spline  fit;  the  third  cubic  passes  through  r = rg  at  F.  S.  360  with  the  same 
slope  as  the  second  cubic  and  through  r = 1,  5 at  F.  S.  430  with  an  assumed  zero 
curvature.  The  radius  distribution  curve  as  composed  of  these  cubics  are  then 
described  by  conic-section  curves  for  input  into  the  computer  programs. 

While  the  lower  blending  radius  could  be  varied,  its  range  would  be  very 
much  limited  in  the  case  of  the  F-4,  Therefore,  a fixed  lower  blending  is 
used  to  assure  a smooth  body  for  the  3DMoC  calculations.  Tids  blending  is 
specified  by  assigning  a fixed  distance  between  points  5 and  F (see  Figure  1). 


4.  FLOW  FIELD  CALCULATION  AND  WAVE  DRAG  REDUCTION 

Once  the  Reometric  variables  and  their  ranges  of  variation  have  been  selected, 
the  fuselage-wing  configurations  corresponding  to  the  nine  cells  of  the  3x3  Latin 
Square  can  be  described.  The  Three-Dimensicxial  Method  of  Chiiracteristics  to- 
gether with  a blunt  body  progi-am  for  providing  the  initial  data  surface  can  then 
be  used  to  calculate  the  flow  fields  around  and  hence  wave  drag  coefficients  of 
these  variation  configurations.  The  wave  drag  coefficients,  in  turn,  can  be  used 
to  determine  the  coefficients  of  the  wave  drag  equations  (1),  from  which  the 
deiiendence  of  the  wave  drag  on  various  geometric  variables  and  the  minimum  wave 
drag  body  can  be  obtained. 

a.  Calculation  of  Flow  Fields 

The  Three-Dimensional  Method  of  Characteristics  program  together  with 
a blunt  body  and  axisymmetric  characteristic  program  was  used  to  calculate  the 
flow  fields  around  the  baseline  F-4  type  fuselage-canopy-wing  configuration  and  the 
variation  configurations  corresponding  to  the  cells  of  the  Latin  Square,  The  wave 
drag  coefficients  were  computed  as  part  of  the  results  to  be  printed  out.  Since  the 
fuselage  nose  is  slightly  blunted,  the  blunt  body  and  axisymmetric  characteristic 
program  was  used  to  provide  a completely  siqjersonic  initial  data  surface  for  the 
Three-Dimensional  Method  of  Characteristics  program  to  proceed.  The  flow-field 
calculations  were  made  at  Mach  2, 5 and  zero  angle  of  atJ  jck.  Since  the  wing  leading 
edge  section  is  described  by  an  ellipse  with  a large  major-to-minor  axes  ratio 
(Equation  14)  the  leading  edge  is  theoretically  blunt.  In  order  to  provide  completely 
supersonic  flow  for  the  characteristics  method  to  calculate,  the  leading  edge  must 
be  subsonic.  At  Mach  2. 5,  a subsonic  leading  edge  has  a sweep  angle  greater  than 
66. 42°.  Hence  a configuration  with  a leading  edge  sweep  of  68°  was  chosen,  which 
is  greater  than  F-4’s  leading  edge  sweep  of  about  51°.  Fortunately,  the  main  geo- 
metric change  that  is  expected  to  yield  appreciable  drag  reduction  is  the  wing-body 


blending  which  does  not  extend  to  the  region  near  the  leading  edge.  Therefore, 
modified  wings  with  considerable  incre:ises  in  the  wing  sweep  can  be  used  for  the 
present  wave  drag  reduction  study  without  ai>preciably  affecting  the  results. 

As  an  example  of  the  flow  fields  calculated  by  using  the  Three-Dimensional 
Method  of  Characteristics,  the  flow  over  the  F-4  type  baseline  at  Mach  2.5  atid  zero 
angle  of  attack  is  presented  here.  Figure  4 shows  the  top  view  of  the  configura- 
tion, the  upper  and  lower  centerline  pressure  distributions  and  the  pressure  along 
the  wing  leading  edge.  Along  the  upper  centerline  of  the  configuration,  the  pressure 
rises  sharply  as  the  flow  hits  the  canopy-fuselage  juncture.  This  signifies  the 
existence  of  an  embedded  shock  wave  created  by  the  canopy.  As  the  flow  spreads 
over  the  canopy  flat  section,  the  pressure  drops  although  the  canopy  profile  is 
almost  strai^t.  A further  drop  of  pressure  is  experienced  as  the  canopy  profile 
curves  back  after  the  flat  section.  It  recovers  to  near  free-stream  pressure  as  the 
w ing  leading  edge  pressure  rises.  Along  the  lower  centerline  the  first  drops 
to  a negative  value  due  to  further  expansion.  It  recovers  somewhat  and  eventually 
comes  close  to  zero,  consistent  with  the  condition  of  zero  angle  of  attack.  The 
wing  begins  with  a high  sweep,  which  decreases  to  a constant  value  of  68*^  near  F.  S. 
190.  Correspondingly,  the  slightly  higher  than  free-stream  pressure  is  obsen'ed  at 
the  beginning  of  the  wing.  Because  of  the  wing-body  interaction,  almost  immediately 
the  wing  leading  edge  pressure  drops  and  does  not  recover  until  after  F.S.  200  when 
the  wing  sweep  drops  to  68°  and  the  larger  wing  span  lessens  the  interaction.  A 
cross  section  of  the  baseline  configuration  at  F.  S.  430  is  shown  in  Figure  5;  the 
wing  has  become  very  thin  aiid  it  is  thicker  hear  the  l^lng  edge  than  near  the  root. 
The  front  view  of  the  baseline  configuration  is  shown  in  Figure  6. 

Of  all  the  variation  configurations,  configuration  6 has  the  maximum  wing-body 
blending.  The  front  view  of  this  configuration  is  shown  in  Figure  7.  A crosf.- 
sectional  view  at  F.S.  350  is  shown  in  Figure  8 together  with  surface  pressure 
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distribution.  At  this  iusohi^o  station,  the  blendinp  is  nearly  maximum,  as  compare<i 
willi  the  baseliiK  cross  scx;tion  shown  by  a dashed  line.  At  F.  j.  ;{.^o  the  blending 
radius  is  still  increasinf’;  1k)wcvoi  , the  rate  of  increase  has  dropped  to  a very 
small  value.  That  menjis  e\i)ansion  has  set  in  upstream  of  F.S.  .'J.'SO,  resultin^j 
in  a neftativo  in  the  blending  region  as  showm. 

b.  Wave  Drag  Equation 

As  shown  in  Equation  (1),  for  a given  Mach  number  and  angle  of  attack,  the  wave 

drag  coefficient  is  assumed  to  be  a quadric  function  of  the  four  chosen  geometric 

variables  x , with  the  coefficients  a.  and  a.,  to  be  determined.  The 

1 4 1 1 j 

reduced  variables  z.,  which  take  the  levels  of  0,  and  * 1 according  to  the  Latin  Squ  ire 
arr.ongement  shown  in  Table  1,  assign  corresponding  values  to  the  geometric 
variables  .x..  For  instance,  according  to  the  first  cell,  z^  and  arc  assigned  level  -1, 
which  corresponds  to  the  minimum  values  of  the  geometric  variables  and  x^  while 
Z2  is  assigned  level  1,  which  corresponds  to  the  maximum  value  of  X2»  and  is 
assigned  level  0 corresponding  to  the  mean  value  of  the  geometric  variable  x^.  Thus 
the  first  cell  specifies  a set  of  values  for  the  geometric  variables  which  in  turn 
defines  a configuration  whose  wave  drag  coefficient  can  then  be  calculated  by  the 
method  of  characteristics  progi*am.  In  this  way,  each  cell  leads  to  one  equation  for 
the  determination  of  the  coefficients  of  the  wave  drag  equation. 

The  wave  drag  coefficients  of  the  baseline  and  the  9 variation  configurations 
are  tabulated  in  Table  3.  These  will  be  used  to  determme  the  wave  di-ag  equation 
in  terms  of  the  reduced  variables  in  the  form 

4 

■^DW  * '=0  ' S ‘ "aV:, 
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It  should  be  noted  that  in  the  original  Latin  Square  surface  fitting  meth'xl 
presented  in  Ileferaice  2,  two  of  the  geometric  variables  have  only  linear  terms 


resulting  in  only  8 coefficients  b.  and  b.^  for  10  equations  and  these  equations  are 

solved  by  a least  square  procedure.  However,  while  the  requirement  of  2 linear 

variables  can  be  relaxed  (see  Appendix  A)  in  the  case  of  5x5  Latin  Sqi  xres,  this 

requirement  does  not  seem  to  apply  to  the  case  of  3 x 3 Latin  Squares.  In  fact,  in 

the  present  case,  these  10  equations  were  used  to  eolve  for  the  set  of  10  unknowns 

•>  2 

b.  and  b.,  that  included  and  z,  in  addition  to  the  origiral  8 coefficients.  The 

i ij  14 

resulting  coefficients  of  the  wave  drag  equation  are  tabula*'ed  in  Table  4. 

c.  Variation  of  Wave  Drag  with  Geometric  Variables 

It  is  instructive  as  well  as  useful  to  represent  the  wave  drag  equation 
graphically.  However,  since  the  wave  drag  depends  on  four  geometric  variables, 
it  is  only  possible  to  show  the  variation  of  the  wave  drag  coefficient  with  respect 
to  two  of  the  variables  while  keeping  the  other  two  variables  constant,  for  instance, 
at  the  baseline  values.  Such  graphs  give  some  ’’feel”  of  the  wave  drag  equation  and 
may  offer  some  insight  about  the  wave  drag  reduction  problem.  Figures  9 to  20 
depict  the  nature  of  the  wave  drag  equation  of  ten  terms  with  coefficients  given  in 
Tabic  4.  These  graphs  were  plotted  on  a Tektronix  equipment.  Each  graph  has 
five  curves  that  correspond  to  five  values  of  the  geometric  variable  shown  at  the 
top  of  the  graph:  thrae  values  equally  divide  the  range  which  is  shown  in  Table  2. 

Figures  9 to  1 1 show  the  variation  of  the  wave  drag  coefficient  with  respect 
to  the  width  a,  which  is  normalized  with  respect  to  the  baseline  width.  As  might 
be  expected,  the  wave  drag  coefficient  increases  with  the  width  of  the  configuration. 
The  curves  are  fairly  straight,  indicating  that  the  wave  drag  is  nearly  proportional 
to  the  width  of  the  configuration,  other  conditions  being  equal.  In  general, 

1'^  increase  of  the  width  increases  the  wave  drag  by  0. 3.5'r,  The  dependence  of  the 

I!) 
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wave  drag  coefficiait  on  the  other  variables  can  be  seen  to  be  far  from  linear 
since  the  distances  between  the  curves  are  quite  different  from  one  another.  The 
dependence  of  the  wave  drag  coefficient  on  the  lower  deck  hieght  h is  shown  in 
Figures  12  to  14.  Since  the  lower  profile  rises  with  increasing  values  of  h,  the 
wave  drag  coefficient  decreases  with  h.  It  is  interesting  to  note  that  the  wave  drag 
attains  a minimum  at  about  h 1.2,  which  represents  a slight  rise  of  the  lower 
profile  from  h = 0.  All  the  curves  are  nearly  flat  near  h = 0,  suggesting  that  a 
slight  change  of  the  position  of  the  lower  profile  has  little  effect  on  the  wave  drag. 

Note  that  the  curves  in  any  of  the  figures  from  Figures  9 to  14  are  parallel  to 
one  another,  because  the  10-term  wave  drag  equation  emtains  no  cross  terms 
of  these  variables  except  r^r.,,  which  represents  the  interaction  between  the  two 
blending  radii  r^^  and  r^.  When  more  configurations  are  included  during  the 
progress  of  the  optimization  procedure,  other  cross  terms  could  be  included, 
as  explained  in  Appendix  A.  However,  the  interaction  between  these  other  vari- 
ables are  expected  to  be  small.  This  observation  was  arrived  heuristically  but 
has  been  verified  by  the  ability  oi  Equation  28  to  accurately  predict  the  minimum 
wave  drag  body  as  will  be  shown.  The  nact  three  figures  15  to  17  show  the  varia- 
tion of  the  wave  drag  coefficient  with  r^,  the  blending  radius  at  F.  S,  280,  for 
different  values  of  one  of  the  other  three  variables  while  the  remaining  two  take  the 
baseline  values.  Similarly,  the  variation  of  the  wave  drag  coefficient  with  Vg,  the 
other  blending  radius  at  F.  S,  3G0,  are  shown  in  Figures  18  to  20.  Although  all 
these  figures  show  a general  increase  of  the  wave  drag  coefficient  with  the  blend- 
ing radius,  there  is  a definite  trend  for  the  cuivcs  to  attain  a minimum  within  the 
ranges.  Most  minima  occur  at  the  lower  side  of  the  radius  scale:  sometimes  it  m;iy 
even  reach  15  inches  as  shown  in  Figure  20.  The  occurence  of  these  minima  is 
significant  because  this  shows  that  the  volume  of  a wing-body  couliguration  can  bo 
increased  by  using  wing-l>ody  blending  without  increasing  the  wave  drag  and  that 
when  the  blending  is  done  propeidy,  the  volume  c:ui  be  increased  with  :m  accompany- 

20 


ing  f eduction  in  the  wave  drag.  Some  interaction  between  and  r,  is  evident  in 
Figures  17  ;md  20,  as  indicated  by  the  cross  over  of  two  or  more  curves. 

d.  Optimization  and  Minimum  Wave  Drag  Body 

The  wave  drag  equation  can  now  be  investigated  to  reduce  the  wave  drag. 

In  the  process  of  determining  the  minimum  wave  drag  configuration,  certain 
geometric  constraints,  such  as  minimum  fuselage  width  or  a given  lower  deck 
height,  must  be  satisfied.  For  each  set  of  geometric  constraints,  there  exists  a 
minimum  wave  drag  configuration.  In  this  section,  the  optimization  procedure  and 
some  minimum  wave  drag  configurations  are  presented  and  discussed.  During 
phase  I study,  a technique  was  developed  (see  Appendix  A)  that  greatly  improves 
the  original  Latin  Square  surface  fitting  method  of  Reference  2,  especially  for 
5 X 5 or  larger  Latin  Squares.  The  improved  optimiz.ation  prcxjcdure  using  the 
improved  Latin  Square  surface  fitting  method  is  verified  in  Appendix  E. 

The  simplest  and  surest  way  to  find  the  minimum  wave  drag  configuration 
subject  to  a given  set  of  geometric  constraints  is  to  use  Equation  13  to  calculate 
the  wave  drag  coefficl«  for  all  allowable  sets  of  levels  of  the  variables  and  pick 
the  set  that  gives  the  least  wave  drag.  The  optimization  procedure  consists  of 
the  following  steps.  * 

1.  Make  a numerical  search  th  rough  the  ranges  of  all  geometric  variables, 
using  the  wave  drag  equation  to  calculate  the  wave  drag  coefficiaits  for 
those  sets  of  levels  that  satisfy  the  constraints. 

2.  Identify  the  set  of  levels  o»  the  variables  that  yields  the  least  wave  drag. 

3.  Prepare  body  description  i-put  data  for  the  minimum  wave  drag  configura- 
tion. 

4.  Use  the  Three-Dimensional  Method  of  Characteristics  program  to  verify 

*The  Mumcrical  Search  Procedure  for  .'ix.'l  Latin  Squares  is  presented  in  Apj)cndix  C: 
it  holds  true  for  3x.3  Latin  Squares  when  the  .space  dimension  is  lowercxl  from  6 to  4. 
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the  prediction  of  the  wave  drag  ecjuation. 


5.  If  the  difference  between  the  predicted  and  calculated  wave  drag  coef- 
ficients exceeds  a certain  criterion,  the  calculated  wave  drag  coefficient 
])rovddes  an  additional  ecjuation  for  the  least  square  fit,  and  the  optimi- 
zation procedure  is  repeated. 

For  the  wave  drag  reduction  problem  of  the  fuselage-wing  configuration, 
as  we  are  concerned  with  in  phase  II,  two  types  of  constraints  are  considered.  The 
first  corresponds  to  setting  one  or  more  of  the  variables  to  a desired  value.  The 
second  corresponds  to  assigning  a fixed  volume  to  the  configuration,  for  instance, 
a certain  percentage  of  the  baseline  volume.  When  two  of  the  variables  are  set 
to  the  baseline  values,  any  of  the  figures  from  Figures  12  to  20  provides  one  or 
more  minimum  wave  drag  bodies.  In  this  respect  these  figures  can  be  quite  useful . 
When  no  constraints  other  than  the  range  limitations  are  imposed,  the  wave  drag 
equation  predicts  a minimum  wave  drag  body  that  corresponds  to  90'7  of  the  base- 
line width,  a raise  of  the  lower  profile  by  1.2  inches  and  a blending  radius  of  3.52 
inches  at  F.S.  280  and  7.63  inches  at  F.  S.  360  (Figure  21).  With  this  configuration, 
a reduction  of  wave  drag  by  4.34^  is  predicted.  The  body  description  input  data 
for  this  minimum  wave  drag  configuration  was  then  prepared  for  the  verification 
run  by  the  Three-Dimensional  method  of  Characteristics  program.  The  results 
showed  a 4. 35'7  reduction  of  wave  drag.  The  difference  between  the  predicted  and 
calculated  wave  drag  reductions  is  within  the  accuracy  of  the  procedure;  therefore, 
the  validity  of  the  10-term  wave  drag  etjuation  is  established.  Figure  22  shows  the 
minimum  wave  drag  configurations  for  various  volume  constraints.  The  percent 
of  wave  drag  increase  was  plotted  against  increasing  volumes  expressed  as  i>crccnt 
of  the  baseline  volume.  The  values  of  the  geometric  variables  that  define  the 
minimum  wave  drag  configurations  are  also  plotted.  It  is  seen  that  a certain  amount 
of  wing-lx>dy  blending  is  present  for  ail  minimum  wave  drag  configurations.  For 
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a given  volume,  the  forward  fuselage  becomes  more  sk/ndcr  to  yield  lower  wave 
drag,  and  the  lost  volume  is  compensated  for  by  the  volume  gained  through  wing-body 
blending. 

e.  Summary  of  Wave  Drag  Reduction  Procedure 

The  present  procedure  can  be  used  to  improve  an  existing  aircraft  or  to  aid 
in  the  design  of  a new  one.  In  the  case  of  improvement,  the  existing  aircraft 
naturally  serves  as  the  baseline.  In  the  case  of  new  design,  the  preliminary  con- 
figuration, which  is  usually  obtained  through  considerations  other  than  the  wave 
drag,  can  serve  as  the  baseline  for  wave  drag  reduction. 

The  next  step  is  to  select  either  the  5x3  Latin  Square  for  the  forward  fuselage 
or  the  3x3  Latin  Square  for  the  blended  wing  configuration  and  to  describe  the  base- 
line and  variation  configurations  using  the  body  description  method  presented  in 
Appendix  C.  The  sections  on  the  description  and  variation  of  configuration  in  Appendix 
A or  the  main  text  should  be  consulted  in  producing  the  body  description.  This  step 
is  time-consuming  but  must  be  done  carefully  to  assure  success  in  wave  drag  cal- 
culations by  the  method  of  characteristics. 

Input  cards  are  then  prepared  according  to  the  instructions  given  in  part  1 and 
2 of  the  user's  manual  (volume  II  of  this  report)  and  fed  to  the  Initial  Value  Surface 
Program  and  the  Three-Dimensional  Metlwd  of  Characteristics  Program  for  wave 
drag  calculations.  Care  needs  to  be  taken  in  the  preparation  of  these  cards,  for  if 
the  calculation  fails  to  proceed  further  the  first  item  to  check  is  the  correctness  of 
the  input  cards. 

After  the  wave  drag  coefficients  have  been  calculated  by  the  Three-Dimensional 
Method  of  Characteristics  Program,  the  Surface  Fit  and  Minimum  Search  Program 
can  be  used  to  define  the  least  wave  drag  configuration  for  a given  set  of  constraints 
as  explained  in  part  3 of  the  user's  manual.  The  program  can  bo  used  as  one  of 
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the  steps  of  a design  procedure  by  providing  the  minimum  wave  drag  body  corre- 
sponding to  const  , lints  imposed  by  other  considerations.  Or  the  program  can  be 
applied  to  generate  a set  of  charts  for  predicting  minimum  wave  drag  Ix)dies 
subject  to  specific'd  constraints.  It  should  be  noted  that  these  charts  are  valid  in 
some  ranges  of  geometric  variables  near  the  baseline.  The  entire  procedure  needs 
to  be  redone  for  a different  baseline  configuration. 
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6.  CONCLUSIONS 


With  regard  to  the  wave  drag  reduction  program,  the  following  remarks 
and  conclusions  can  be  made. 

a.  The  present  optimization  procedure  developed  in  the  phase  I study  is  useful 
and  versatile.  It  can  be  used  for  other  optimization  purposes. 

b.  Together  with  the  Three-Dimaisional  Method  of  Characteristics,  the  procedure 
can  be  used  to  obtain  the  minimum  wave  drag  configuration  for  designing  new 
airplanes  or  modifying  existing  ones. 

c.  The  optimization  procedure  has  been  verified  through  applications  to  the  F-4 
forward  fuselage  using  a 5 x 5 Latin  Square  and  to  an  F-4  type  fuselage-wing 
configuration  using  a 3 x 3 Latin  Square.  The  procedure  has  also  proven  Itself 
by  correctly  predicting  the  von  Karman  ogive  as  the  minimum  wave  drag  body 
for  a given  configuration  length  and  base. 

d.  Within  the  ranges  of  variation  of  the  chosen  geometric  variables,  the  following 
rule-of-thumb  precentage  reductions  of  the  wave  drag  are  obtained.  In  the 
case  of  the  F-4  fuselage,  forj  every  inch  the  nose  is  lengthened,  the  wave  drag 
is  reduced  by  sll^tly  over  one  percent,  and  for  every  percent  the  fuselage 
volume  Is  decreased,  the  wave  drag  is  reduced  by  about  three  quarters  of  a per- 
cent. In  the  case  of  the  F-4  type  fuselsage-wing  configuration,  the  wave  drag  is 
reduced  by  about  half  a percent  for  each  percent  the  blended-wlng  configuration 
is  narrowed. 

e.  The  present  application  to  wave  drag  reduction  is  limited  only  in  the  capability 

i 

of  the  wave  drag  computational  techniques.  First,  a completely  supersonic  flow 

field  is  required  for  the  characteristics  method  to  be  applicable.  For  a given 

' ] 

configuration,  this  requirement  sets  a lower  limit  on  the  free-stream  Mach 
number.  Secondly,  the  computer  program  at  the  present  stage  of  development 

i 

requires  that  all  corners  and  ^ges  be  faired  with  smooth  curves  to  yield  a 
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unique  surface  normal  everywhere  on  the  configuration. 


At  free-stream  Mach  niunbers  below  the  lower  limit  stated  in  e. , subsonic 
regions  would  occur  at  the  fuselage-canopy  juncture  of  the  configuration  under 
consideration.  Further  studies  on  calculations  of  local  subsonic  regions  are 
needed  to  provide  methods  for  supplementing  to-day's  wave  drag  computational 
techniques  at  lower  free-stream  Mach  numbers. 


26 


REFERENCES 


1.  Chu,  C.  W. , Der,  J.  Jr.,  and  Ziegler,  H. , "Wave  Drag  Reduction  for  Aircraft 
Fuselages,"  Interim  Report  NOR  75-70,  August  1975,  Northrop  Corporation, 
Hawthorne,  California, 

2.  Redlich,  O. , and  Watson,  F,  R. , "On  Programs  for  Tests  Involving  Several 
Variables,"  Aeronautical  Engineering  Review,  Vol,  12,  No,  6,  June  1953, 
pp.  51-59. 

3.  Chu,  C.  W. , "Compatibility  Relations  and  a Generalized  Finite-Difference 
Approximation  for  Three-Dimensional  Steady  Supersonic  Flow,"  AIAA  Journal 
Vol.  5,  No.  3,  March  1967,  pp.  493-501. 

4.  Chu,  C.  W.  and  Powers,  S.  A. , 'The  Calculation  of  Three-Dimensional 
Supersonic  Flows  Around  Spherically-Capped  Smooth  Bodies  and  Wings," 
AFFDL-TR-72-91,  Vol.  I,  Theory  and  Applications,  September  1972. 

5.  Chu,  C.  W. , "Calculation  of  Three-Dimensional  Supersonic  Flow  Fields 
about  Aircraft  Fuselages  and  Wings  at  General  Angies  of  Attack,"  NOR  72- 
182,  March  1973,  Northrop  Corporation,  Hawthorne,  California. 

6.  Chu,  C.  W. , and  Powers,  S.  A.,  "Determination  of  Space  Shuttle  Flow  Field 
by  the  Three-Dimensional  Method  of  Characteristics, " TMX  2506,  Feb.  1972, 
pp.  47-63,  NASA. 

7.  Chu,  C.  W. , "A  New  Algorithm  for  Three-Dimensional  Method  of  Character- 
istics," AIAA  Journal,  Vol.  10,  No.  11,  November  1972,  pp.  1548-1550. 

8.  Chu,  C.  W. , "Calculation  of  Supersonic  Flow  Fields  about  Slab  Delta  Wings 
and  Space  Shuttle  Wing-Body  Ccmfigurations, " NOR  73-007,  April  1973, 

Northrop  Corporation,  Hawthorne,  California. 

9.  Chu,  C.  W. , "Supersonic  Flow  About  Slab  Delta  Wings  and  Wing-Body  Configura- 
tions," Journal  of  Spacecraft  and  Rockets,  Vol.  10,  Na  11,  November  1973, 

pp.  741-742. 


27 


Best  Available  Copy 

TABLE  1.  Latin  Square  Arrangement 


z 

2 

3 

CELL  Code: 


1 0 

®n 

0 -1 

-1  1 

0 -1 

® . 

-1  0 

-1  -1 

© 

0 0 
0 0 

1 0 
1 1 

0 i 

© 

0 1 
1 -1 

© 

1 

o 

TABLE  2.  Variables,  Ranges,  and  Baseline  Values 


Reduced 

Variables 

Geometric 

Variables 

Ranges  of  x. 

[ Baseline  Values 

• 

mmm 

WmSm 

X. 

max 

z. 

1 

X. 

i 

a 

0.9 

1.15 

-0.2 

1.0 

*2 

^1 

1.5 

36. 30 

"3 

^2 

1.5 

61.30 

-1.0 

^4 

h 

3.0 

-6.0 

-0,333 

0.0 
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TABLE  3.  Wave  Drag  Coefficients 
(Base  on  Wing  Area  of  530  Sq.  Ft.) 


0.0104602 

0.0116525 

0.0107868 

0.0094555 

0.0104013 

0.0128320 

0.0116629 

0.0114297 

0.0114309 

Baseline  Value  = 0.00979359 


TABLE  4.  Coefficients  of  Wave  Drag  Equation 


b =0.0104013 
o 

Values  of  b.  x 10^ 


i 

1 

2 

3 

4 

bj  X 10^ 

5. 7851670 

3. 638333 

7.459000 

2. 706667 

Values  of  b.^  x 10^ 


1 

2 

3 

4 

1 

0.6171822 

2 

2. 598000 

-1.181682 

3 

5.391000 

4 

2.227318 
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rv-  «;ectlons  Description 

Fuselage  and  Wing  Cross  Sectio 
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FIGURE  2.  Geometric  Variables  for  Fuselage  Variation  and  Wing-Body  Blending 
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FIGURE  5.  Cross  Section  ol  F-4  Ty\yo  Baseline  at  Fuselage  Station  430 


FIGURE  6.  Front  VicTvs  of  F-4  Type  Baseline  Cross  Sections 
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FIGUHE  7.  Front  Views  «»f  Confi};ur:Uion  C Cross  Sections 


FUSELAGE  CROSS  SECTION 


37 


FIGURE  8.  Cross  Section  and  Pressure  Distribution  at  F.  S.  350  of  Blended  Wing  Configuration  6 
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FIGURE  11.  Variation  of  Wave  l>raR  Coefficient 
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FIGURE  14.  Variation  of  Wave  Drag  Coefficient  ^ ~ ~ 
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FIGHUE  If).  Variation  ot  Wave  Orat;  Co«*Hieient  vs 
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FIGURE  16.  Variation  of  Wave  Drag  Coefficient  vs 
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FKUIUE  17.  Vai  iati«m  ul  Wave  Drat;  Coeflieicnl  vs 
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FlGt'liE  20.  Variation  of  Wave  Drag  Coefficient  vs  i'2  ^ 
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FIGURE  22.  Minimum  Wave  Drag  Ccnfiguratlons  for 
I Various  Volumes  of  Baseline  Volume) 
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APPENDIX  A 


WAVS  DRAG  REDUCTION  FOR  FORWAP">  FUSELAGES 
1.  INTRODUCTION 

The  complete  wave  drag  reduction  program  has  beas  carried  out  in  two  phases. 
It  the  first  phase,  a procedure  was  developed  for  minimizing  the  wave  drag  of  a 
forward  fuselage  and  canopy  configuration.  In  the  second  phase,  the  procedure  was 
expanded  to  account  for  the  influence  of  the  wing  and  wing-body  blending  on  the  over- 
all wave  drag.  In  this  appendix,  the  research  performed  in  the  first  phase  of  the 
program  is  presented.  The  basic  approach  is  given  first,  which  consists  of  the  for- 
mulation of  the  problem  and  a brief  account  of  the  two  basic  methods.  It  is  followed 

by  a discussion  of  the  method  of  describing  the  body  and  the  selection  of  geometric 
variables  and  their  ranges  for  defining  a family  of  configurations.  Then  the  flow  field 

calculation  and  the  wave  drag  equation  are  presented.  Sample  results  of  the  calculated 
flow  fields  are  given,  and  calculated  wave  drag  coefficients  are  tabulated.  These 
coefficients  were  used  to  derive  the  wave  drag  equation  which  expresses  the  wave 
drag  as  a function  of  the  geometric  variables.  A new  concept  which  enabl.;5  the  wave 

drag  equation  to  "learn”  from  experience  to  improve  its  performance  is  also  pre- 
sented. This  concept  proved  very  useful  in  achieving  successftil  results  during  the 

Phase  I work  and  can  be  applied  to  improve  other  optimization  procedures  using  Latin 
Square  sampling.  Then  the  wave  drag  reduction  procedure  and  the  types  of  geometric 
constraints  imposed  by  design  requirements  considered  in  the  procedure  are  pre- 
sented, followed  by  the  discussion  of  the  wave  drag  equation  and  the  new  concept  for 

improvement.  Various  aspects  of  the  wave  drag  reduction  procedure  are  demon- 
strated using  the  F-4  fuselage  aa  the  baseline;  the  results  are  presented  and  dis- 
cussed. Some  characteristics  of  the  wave  drag  equation  are  plotted,  and  concluding 
remarks  are  made  for  the  Phase  I work. 
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■2.  A1’P1U)AC!1 


In  this  section  the  formulation  of  the  wave  drag  reduction  problem  is  outlined, 
and  the  two  basic  methods  to  be  used  in  this  study  ai-e  introduced. 


a.  Formulation 

'I'he  baseline  configuration  can  bo  described  by  a set  of  geometric  variables.  A 
family  of  configurations  including  the  baseline  can  be  generated  by  assigning  different 
values  to  some  or  ail  of  the  geometric  variables.  If  the  wave  drag  can  be  expressed 
as  a function  of  these  variables,  a particular  set  of  values  of  these  variables  that  gives 
the  least  wave  drag  can  be  found  by  minimizing  the  function.  This  set  of  values  then 
produces  the  minimum  wave  drag  configuration. 

The  key  to  this  problem  is  how  to  obtain  such  a functional  expression  for  the 
wave  drag.  In  this  study,  sLx  geometric  variables  are  considered  in  defining  the  family 
of  configurations.  If  each  variable  assumes  five  values,  the  evaluation  of  the  partial 
derivatives  with  respect  to  these  six  variables  for  a Taylor-series  type  expression 
would  require  5*’  or  wave  drag  calculations.  This  is  obviously  not  feasible.  In 

the  present  approach  the  wave  drug  coefficient  is  first  assumed  to  Ije  of  the  form 

G 5 5 


C„.„  =»  a + V'  a.x.  ■* 

I)W  o I I 

1=^1  j=i  i=2 


V y*  a..x.x. 


(1) 


where  x.  (i  = 1,0)  arc  the  geometric  variables  and  a^  and  a.,  are  to  be  determined. 

The  Latin  Square  sampling  technique^  is  used  to  sample  25  sets  of  values  of  x^  out  of 
the  total  population  of  15,  025  sots.  The  wave  drag  coefficients  for  the  configurations 
defined  by  each  of  the  25  sets  are  then  calculated  by  the  Three-*  Dimensional  Metliod  of 


Characteristics. 


2,5. 


When  the  25  calculxitcd  wave  drag  coefficients  and  the 


corres  nding  values  of  the  geometric  variables  Xj  arc  substituted  into  Eq.  (1),  25 
linear  c-quations  for  the  17  unknowns  a.  and  a.^  arc  obtained.  A least-squares  procedure 
is  used  to  solve  these  equations  for  the  a.  and  a.^,  which  arc  then  substituted  back  into 
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Eq.  (1)  to  produce  the  functional  expression  for  the  wave  drag  in  terms  of  the 
geometric  variables.  By  minimizing  in  Eq.  (1)  subject  to  a given  set  of  con- 
straints, e.  g. , a given  length  and  width  of  the  fuselage,  the  minimum  wave  drag  con- 
figuration corresponding  to  the  given  set  of  constraints  is  determined.  In  the  present 
study,  a numerical  search  pro<;edure  is  used  to  find  the  minimum  wave  drag  configura- 
tion (Appendix  D). 

b.  Latin  Square  Sampling  Technique 

The  Latin  Square  method,  which  has  mostly  been  used  in  agriculture  and 
Hological  research,  is  a very  efficient  sampling  technique  and  is  much  better  than 
random  sampling.  For  this  study  a particular  type  of  Latin  Square  (the  orthogonal 
squares)  suitable  for  a variety  of  technical  problems^  is  adopted.  With  this  type  of 
Latin  Square  arrangement,  a 5X5  square  is  the  correct  size  for  six  geometric 
variables  Xj,  each  taking  five  values.  It  is  convenient  to  introduce  the  reduced 
variables  z.,  which  are  related  to  the  geometric  variable  through 

X,  -t-  X.  ~ \ 

max  ^min  * \ ^max  ^min/  19  \ 

= ^ (2) 

where  subscripts  max  and  min  denote  the  maximum  and  minimum  values,  respectively. 

Corresponding  to  the  five  values  of  the  geometric  variables  x^,  the  reduced  variables 

Zj  always  assume  the  five  levels  of  0,  *1,  i2.  The  5X5  Latin  Square  arrangement  in 

terms  of  the  levels  of  the  reduced  variables  is  shown  in  Table  1.*  It  is  seen  that  in  this 

way  the  Latin  Square  arrangement  remains  the  same  whatever  the  values  of  the 

geometric  variables.  At  first  glance,  it  appears  that  the  roles  of  z,  and  z_  are 

i o 

unique  since  their  levels  are  arranged  regulariy  in  Table  1 and  their  nonlinear  terms 
are  excluded  from  Eq.  (1).  This  is  true  for  the  conventional  approach.  However,  it  is 
shown  in  Appendix  B i.hat  by  rearranging  Table  I any  pair  of  the  reduced  variables  can 
have  their  levels  arranged  regularly  as  Zj  and  7.^.,  and  it  will  be  shown  later  in  this 

• A discussion  of  the  Latin  Square  construction  is  given  in  Appendix  B. 
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appendix  that  some  or  all  nonlinear  terms  of  Zj  and  Zg  can  be  included  in  the  wave  drag 
equation. 

c.  Three-Dimensional  Method  of  Characteristics 

2 3 

In  this  study  the  Three-Dimensional  Method  of  Characteristics  * is  used  to 
calculate  the  wave  drag  coefficients  for  the  fuselage-canopy  configurations  sampled 
by  the  Latin  Square  technique. 


/ 
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3.  DKSCRIPTION  AND  VARIATION'  OF  FORWARD  FUSKlj\GE 

Three-dimensional  body  description  requires  a great  deal  of  effort  which  at 
times  becomes  extremely  tedious.  Basically,  two  types  of  description  can  be  made, 
analytical  and  numerical.  The  latter  can  describe  complicated  geometry  accurately, 
but  is  cumbersome  for  preparing  input  data.  The  former  is  much  simpler  to  input  and 
can  define  a family  of  configurations  based  on  a few  geometric  variables.  Hence, 
analytical  description  is  chosen  for  the  present  study. 

a.  Body  Description 

Every  configuration  has  a number  of  generating  lines,  such  as  the  upper  profile, 
the  lower  profile,  the  maximum  breadth  line,  or  the  wing  leading  edge.  In  the  present 
body  description  procedure,  each  generating  line  is  divided  into  a number  of  segments 
to  permit  each  segment  to  be  described  by  a conic-section  curve.  At  each  cross  sec- 
tion of  the  configuration,  simple  analytic  curx'es,  e.g. , the  ellipse  or  cubic,  connect 
any  two  adjacent  generating  lines  to  form  the  contour  of  the  cross  section.  The  con- 
figuration is  thus  described  analytically  by  simple  low-order  curves.  For  a smooth 
body  a unique  normal  to  the  surface  exists  everywhere,  and  this  condition  usually 
requires  slope  continuity  at  the  junctures  between  two  contour  curves  or  two  segments 
of  a generating  line. 

The  fuselage  Is  located  in  a right-handed  coordinate  system  where  the  Y-axis  is 
aligned  with  the  fuselage  axis;  the  X-axis  is  spanwise  and  the  Z-axls  is  up.  A schematic 
of  the  fuselage-canopy  configuration  is  shown  in  Figure  i ano  u typical  cross  section, 
in  Figure  2.  All  generating  lines  are  represented  by  a g.  ;ral  curve  fit  of  conic  sec- 
tions in  several  segments.  The  conic-section  curves  la'f  h • f ..'TI 

( M a P.Y  + Q.  i (r.y^  I (3) 

\x.y  1 1 

where  j » 1 , 5.  A straight  line  is  a special  case  with  R = S * T * 0.  Each  curve  can 
be  divided  into  as  many  segments  as  necessary  to  provide  adequate  body  description. 
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Kacli  segment  must  be  continuous  with  the  previous  segment  and  with  very  few  excep- 
tions the  slope  must  be  continuous  at  the  junctures  to  satisfy  the  requirement  of  a 
unique  normal  to  the  surface.  I’nlilvC  others,  line  ti  is  not  a true  generating  line.  It 
is  actually  a shape  factor  which  expresses  the  distance  to  a L -^ent  line  as  a function  of 

Y and  is  discussed  in  Section  3b.  For  lines  7 and  8,  only  and  are  fitted  by  conic 
section  curves  while  X„  and  X are  given  by  Eq.  (10a)  in  Section  3c. 

< O 

As  shown  in  Figure  2,  the  contour  of  the  cross  section  begins  with  a straight 
line  representing  the  canopy  flat.  The  canopy  contour  from  points  4 to  5 is  circular 
but  can  be  elliptic  in  general.  The  upper  hiselage  is  rei)resented  by  an  elliptic  curve 
from  points  1 to  2.  The  lower  fuselage  is  represented  by  a general  conic- section 
curve  from  points  2 to  3;  the  bottom  Qat,  by  a straight  line  from  point  3 to  the  center- 
line.  The  intersection  between  the  canopy  and  the  fuselage  is  faired  by  a cubic  from 
points  7 to  8.  The  derivation  of  the  cubic  equation  is  presented  in  Section  3c.  The 
equation  for  the  canopy  is  given  by 

2 2 

■ Z - Zg  (Y)  ■ 

Z“  (Y)  - (Y) 

If  it  is  a circular  arc,  then  (Z^  - Z^)^  * (Xg  - X^)^.  The  equation  for  the  upper 
fuselage  is 

2 2 


■ z - Z2  (Y) 

1 

■ X - Xj  (Y) 

Z^  (Y)  - Z2  (Y) 

T 

X2  (Y)  - X^  (Y) 

X - (Y) 

[xg  (Y)  - X^  (Y)J 

-1*0 


(4) 


b.  Shape  Factor 

Referring  to  Figure  2,  the  contour  curve  from  points  2 to  3 has  a zero  slope  at 
point  3 and  the  slope  approaches  infinity  at  point  2.  As  in  aircraft  lofting  practice, 
the  shape  of  this  contour  curve  is  determined  by  specifying  the  distance  from  the  origin 
to  one  of  its  tangent  lines  that  makes  a 45“  angle  with  the  x-axis.  Hence,  this  distance , 
which  is  designated  by  b in  the  following  derivation  of  the  equation  for  the  contour 
curve,  may  be  regarded  as  a shape  factor  for  this  curve. 
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This  equation  for  the  45''  tangent  line  is 

X = Z - /2  b 

which,  when  the  origin  is  translated  to  point  3,  becomes 

X'  = Z’  + H (G) 

where  Z’  = Z-Z^,  X’  = X-X^,  and  H = Z^-X^  ~/2  b.  In  the  new  coordinate  system, 
the  general  quadric  equation,  satisfying  the  conditions  Z'  (X’2)  = Z’^,  X’  (0)  = 0, 

,(^)2  * 

K(Z^  X'  - X^Z’J^  + Z’  (XV-  X^)=  0 (7) 

which  represents  a family  of  curves  with  X as  a parameter.  To  determine  K. 

Eq.  (6)  is  substituted  into  Eq.  (7)  to  yield  a quadratic  equation  in  the  form 

AZ'^  + BZ'  + C « 0 (8) 

and  the  condition  for  Eq.  (6)  to  be  a tangent  line  is  that  Eq.  (8)  should  have  a double 
root;  i.  e. , - 4AC  * 0,  which  leads  to 

(H-X^)2 

^ “aZ^X'2(H+  Z^  - 

Equation  (7)  with  K given  by  Eq.  (9)  is  the  equation  for  the  contour  curve  from 
points  2 to  3.  The  range  of  variation  of  this  curve  obtainable  by  applying  this  equation 
is  illustrated  in  Figure  3,  where  a family  of  conic- section  curves  is  given  for  dif- 
ferent values  of  b.  By  increasing  the  distance  b,  the  curve  is  seen  to  vary  from 
almost  a straight  line  to  a sharply  bending  curve  approaching  the  two  sides  of  a right 
triangle. 

c.  Fairing  Curve 

The  intersection  between  the  canopy  and  fuselage  is  faired  by  a cubic  from  points 
7 to  8 (Figure  2).  The  projections  Z,^(Y)  and  Zg(Y)  of  lines  7 and  8 on  the  Y-Z  plane 
are  given  by  Eq.  (3)  where  the  coefficients  P,  Q,  R,  S,  T are  input  quantities.  The 
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proji'ctioos  X,^  (Y)  and  (Y)  on  the  X-Y  plane  arc  obtained  by  solving  Eqs  (4)  and  (5). 
respectively. 


>C^(Y)=X^M)C-X^) 


Xg  (Y)  = X,  4 (X^  . X,) 


/ Z_  - 7.. 

1 - 7 2. 


1 


7. , - 
•4  a 


Z - Z- 

_8 ^ 

Z.  - Z„ 


1/2 


1/2 


(10a) 


1 


'1  '*2  1 

The  fairing  curve  matches  the  slopes  of  the  ellipses  at  both  end  points  7 and  8.  The 
slopes  are  obtained  by  differentiation  of  Eqs.  (4)  and  (5) 

-X,'" 


■a  • -tefe:)’ 


(lOb) 


X-giY) 


■'8  “8  1 1 
The  cubic  equation  that  satisfies  tlic  first  conditions  of  Eqs.  (10a)  and  (10b)  can  be 
written 


X = X? 


(11) 


The  coefficients  c and  d are  obtained  by  applying  the  last  conditions  of  Eqs.  (10a)  and 
(10b) 

c = .T  (Xg  - X.)  - (X;  + 2X^)  (Zg  - Z^)  (12) 

d=-2(Xg-X^)  ^ (X;,  4 x^)(Zg*Z^) 

Equation  (11)  with  c and  d given  by  Eq.  (12)  is  then  the  cubic  equation  for  the  fairing 
curve  from  points  7 to  8.  'Ihe  quantities  X^,  Xg,  X’^  and  X’g  are  given  by  Eqs.  (Toa) 
and  (10b)  while  Xj,  X.,,  Xj,  X^.,  Zj,  Z^,  Zj,  Zj.,  Z^  .and  Zg  arc  obtainwl  from  Eq. 

(.3)  where  coefficients  1’^,  Q^,  i(,,  and  T are  input  (|uantities. 
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d.  Cicometric  Variables  and  Itan^cs  of  \'ariation 


Six  geometric  variables  were  selected  to  generate  a family  of  fonvard  fuselages 
including  the  F— 1 baseline.  .\s  illustrated  in  Figure  4,  these  variables  are  the  length 
; of  the  fuselage  nose,  the  horizontal  displacement  a and  vertical  displacement  a'  of 

the  maximum  horizontal  breadth  line,  the  shape  factor  b,  the  lower  deck  height,  h, 
and  the  bottom  flat  width  f.  As  discussed  previously  the  shape  factor  b is  the  distance 
from  the  origin  to  the  15'  tangent  line  which  is  tangent  to  a contour  curve  of  the  lower 
fttselagc  cross  section.  All  configurations  generated  by  these  variables  have  the  same 
canopy  as  the  F-4  and  must  satisfy  the  over-lhe-nose  and  over*the-side  view  line 
limitations. 

The  correspondence  between  the  geometric  variables  and  the  reduced  variables 
defined  in  Eq.  (2),  the  ranges  of  variation  of  these  geometric  variables,  and  the  baise- 
line  values  of  these  variables  are  tabulated  in  Table  2.  The  variety  of  configurations 
that  can  be  obtained  using  these  geometric  variables  within  their  ranges  of  variation 
is  illustrated  in  Figures  5 and  0.  Figure  3 shows  the  top  and  side  views  of  the  con- 
figurations. The  solid  lines  show  the  baseline  configuration  of  the  F-4.  The  dash-dot 
lines  describe  two  extreme  variations  of  the  fuselage  shape.  The  variation  of  the  fuse- 
lage shape  ahead  of  the  canopy  has  to  satisfy  the  over- thc-nose  view  constraint,  and 
the  width  of  the  fuselage  is  subject  to  the  over-thc-side  view  constraint.  The  ratlar 
installation  Imiioses  a minimum  requirement  for  the  width  and  limits  the  amount  by 
which  the  lower  profile  can  be  moved  upward.  Notice  that  the  figiirc  shows  only  two 
extreme  examples  of  the  variations.  The  long  nose  version  can  be  combined  with  a 
wide  fuselage  or  the  short  nose  version  can  be  combined  with  a slender  fuselage  to 
yield  other  intermediate  conficurations.  Figure  (>  demonstrates  the  type  of  configura- 
tion variation  that  can  be  achieved  through  the  chosen  variables.  Cross-sectional 
views  at  the  three  hiselage  stations  indicated  in  Figure  5 are  shbwn  iji  Figure  G.  The 
solid  lines  describe  the  baseline  configuration  while  the  dashed  lines  describe  two 
extreme  variation?  of  lliecros.s  sections.  The  circle  at  F..S.  50,  which  represents 


TABLE  2.  VARIABLES,  R^VNCES,  AND  BASELINE  VALUES 


Reduced 

Variables 

"i 

Geome  trie 
Variables 

X. 

1 

Ranges  of  x. 

Baseline  Values 

X. 

^min 

X. 

L 

max 

2 . 
1 

X. 

1 

t 

84 

104 

-1.38 

, 87.1 

^2 

b 

-2/3 

2/3 

0 

* 

"3 

f 

0. 

1.5 

2/3 

1. 

h 

5. 

-10. 

-2/3 

0. 

^5 

a’ 

0.6 

1.2 

2/3 

1. 

a 

0.9 

1.15 

-0.4 

1. 

*Baseline  values  for  x.  vary  from  station  to  station,  x.  is  taken 

rn  i ri 

to  be  the  Local  baseline  value  minus  2/3  of  the  difference  between  the 
baseline  value  and  the  local  minimum  value  and  x.  is  taken  to  be  the 

i 

max 

local  baseline  value  plus  2/3  of  the  difference  between  the  baseline  value 
and  the  local  maximum  value. 
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FIGUHE  5.  VAHlA'liON  OF  FISFLAGF  PROFILES  AND  BREADTH 


ARIATION  OF  FUSEI.AOt;  CROSS  SECTIONS 


the  minimum  area  for  the  radar,  limits  the  amount  by  which  the  maximum  breadth  line 
can  move  in  or  the  lower  profile  can  move  up.  The  over-the-side  view  line  governs 
how  far  the  maximum  breadth  line  can  move  out  and  up.  Again  only  two  extreme 
examples  are  given.  Other  configurations  can  be  easily  constructed.  The  important 
variable  for  the  cross-sectional  variation  is  the  shape  factor  b which  drastically  alters 
the  shape  of  the  lower  fuselage  between  the  maximum  breadth  line  and  the  lower  sur- 
face flat. 

It  should  be  fairly  obvious  that  the  choice  of  geometric  variables  and  their  ranges 
of  variation  are  somewhat  arbitrary.  The  present  selection  allows  a fairly  wide  range 
of  variation  of  configurations  which  are  similar  to  the  baseline  configuration.  Other 
selections  can  be  made  such  that  the  range  of  variation  is  more  limited  or  the  family 
of  configurations  is  less  similar  to  the  baseline. 

Some  guidelines  for  the  selection  of  geometric  variables  and  their  ranges  can 
be  given. 

(1)  The  two  linear  variables  should  be  selected  if  their  effects  on  the  wave  drag 
are  small  or  if  they  affect  the  wave  drag  nearly  linearly  and  their  interac- 
tions with  other  variables  are  very  small.  The  interaction  between  two 
variables  is  the  influence  of  either  variable  on  the  wave  drag  contribution 
of  the  other.  Nonlinear  effects,  however,  can  be  included  through  a new 
concept  which  will  be  discussed  in  Section  4c. 

(2)  The  four  remaining  nonlinear  variables  should  be  selected  in  such  a way 
that  their  interactions  with  each  other  are  small  or  nearly  linear. 

(3)  The  range  of  each  variable  should  be  as  small  as  necessary,  just  wide 
enough  to  serve  the  particular  purpose. 

(4)  In  general  the  nonlinear  variables  and  their  ranges  should  be  selected  in 
such  a way  that  their  influences  on  the  wave  drag  are  of  the  same  order  of 
magnitude. 
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4.  1 LOW  ITKI.D  CALCULATION  AND  WAVK  DHAG  KQUATION 


The  wave  drag  eoetTieionts  of  the  baseline  configuratiem  and  the  25  configurations 
pickc'd  by  the  Latin  Square  teclinique  are  calculated  by  applying  the  Threc-Uimensio»  <i 
Method  of  Characteristics.  I'he  resulting  wave  drag  coefficients  are  then  uf  • ■ lO  ob- 
tain the  wave  drag  equation. 

a.  Mow  Over  I'-  1 Baseline 

As  an  example  of  the  flow  fields  calculated  by  using  the  rhrce-Dimensional 
Method  of  Characteristics,  the  flo..'  over  the  K-1  baseline  at  Mach  2.5  and  zero  angle 
of  attack  is  presented  here,  figure  7 shows  the  top  and  front  views  of  the  configura- 
tion, four  cross-sectional  views,  and  the  upper  and  lower  centerline  pressure  distri- 
butions. Along  the  upper  centerline  of  the  configuration,  the  pressure  rises  sharply 
as  the  flow  hits  the  canopy- fuselage  junctm*c.  This  signifies  the  existence  of  an 
embedded  shock  wave  created  by  the  canopy.  At  the  pressure  peak  the  local  Mach 
number  drops  to  1. 17.  As  the  flow  spreads  over  the  canopy  flat  section,  the  pressure 
drops  although  the  canopy  profile  is  almost  straight.  A sudden  drop  of  pressure  is 
experienced  as  the  canopy  profile  curves  back  after  the  flat  section.  Along  the  lower 
centerline  the  first  drops  to  a negative  value  due  to  further  expansion.  It  recovers 
somewhat  and  eventually  comes  close  to  zero,  consistent  with  the  condition  of  zero 
angle  of  attack. 

Figure  H shows  polar  plots  of  the  pressure  distributions  at  four  fuselage  stations 

with  cross  sections  shown  in  Figure  7.  Since  the  forward  fuselage  ahead  of  the  canopy 

droops  somewhat,  the  pressure  is  higher  on  the  upper  fuselage  than  that  on  the  lower 

fiiselage  at  zero  angle  of  attack.  At  F.S.  :i5.4  the  increases  monotonically  from  a 

slightly  negative  value  at  the  bottom  to  a positive  value  at  the  top.  A noticeable 

pressure  rise  can  be  seen  at  F..S.  85.3  near  the  fuselage-canopy  juncture.  A similar 

pressure  peak  is  visible  at  F.S.  145,  at  which  C has  become  negative  e.xcept  near 

P 
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FUSELAGE  STATION  Y 


FIGURE  7.  F-4  FUSELAGE  AND  CALCULATED  CENTERLINE 
PRESSURE  DISTRIBUTIONS 
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(a)  FS  35.4 


(b)  Fb  85.3 


FIGUUK  8.  POLAR  GRAPHS  OF  PRFSSURK  DISTRIBUTIONS  AT  FOUR 

FUSELAGE  STATIONS  | 
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the  bottom  flat.  In  general  the  flow  continues  to  expand  downstream.  By  F.S.  203  the 
pressure  everywhere  has  become  slightly  less  than  the  free-stream  pressure,  and  the 
variation  of  Cp  is  small  and  nearly  monotonic. 


b.  Calculation  of  Wave  Drag  Coefficients 

As  shown  in  Eq.  (1),  for  a given  Mach  number  and  angle  of  attack,  the  wave 
drag  coefficient  C^^is  assumed  to  be  a quadric  hinction  of  the  six  chosen  geometric 


variables  x^, 


. . X.  with  the  coefficients  a.  and  a,,  to  be  determined.  The  reduced 
b 1 ij 


variables  z..  which  take  the  levels  of  0.  tl,  i2  according  to  the  Latin  Square  arrange- 
ment shown  in  Table  1,  assign  corresponding  values  to  the  geometric  variables  x.. 

For  instance,  according  to  the  first  cell,  and  Zg  are  assigned  level  -2,  which  cor- 
responds to  the  minimum  values  of  the  geometric  variables  Xj  and  Xg  while  z^  is 
assigned  level  2,  which  corresponds  to  the  maximum  value  of  x^.  The  other  three 
reduced  variables  are  assigned  levels  corresponding  to  intermediate  values  of  the 
geometric  variables.  Thus  the  first  cell  specifies  a set  of  values  for  the  geometric 
variables  which  in  turn  defines  a configuration  whose  wave  drag  coefficient  can  then 
be  calculated  by  the  method  of  characteristics  program.  In  this  way,  each  cell  leads 
to  one  equation  for  the  determination  of  the  a.  and  a.^  in  the  wave  drag  equation. 

Since  the  method  of  characteristics  is  based  on  hyperbolic  equaticns,  the  com- 
plete flow  field  must  be  supersonic  in  order  to  apply  the  method.  As  mentioned  above, 
the  local  Mach  number  dropped  to  . 17  near  the  fuselage  canopy  juncture  of  the  F-4 
at  free-stream  Mach  2.5  and  zero  angle  of  attack.  Local  subsonic  regions  would  likely 
occur  for  some  of  the  25  variation  configurations  if  the  free-stream  Mach  number 
should  be  lowered.  Hence  the  wave  drag  coefficients  for  the  F-4  fiiselage  and  the 
25  variations,  each  corresponding  to  a cell  of  the  Latin  Square,  were  calculated  using 
the  Three-Dimensional  Method  of  Characteristics  program  at  Mach  2. 5 and  zero  angle 
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of  attack.  The  calculated  wavi*  drajj  coc f fie ic tits  (based  on  a win^  area  of  5.30  sq.  ft. ) 
are  tabulated  in  Table  The  mean  value  of  the  3.5  coefficients  is  slightly  higher  than 
the  baseline  coefficient. 

c.  Wave  Drag  figuation 

In  the  wave  drag  equation,  Kt(.  (1),  two  of  the  geometric  variable.s,  i.e. , Xj  and 
X,.,  have  onlv  linear  terms  while  other  variaides  have  nonline.'’r  terms.  In  this  wav 
there  are  onlv  17  eoefficients  a.  and  a.,  while  there  are  25  equations  availalile  for  their 

t ij 

evaluation.  It  therefore  seems  logical  to  include  nonlinear  terms  of  x^  and  x^,  in  the 

equation  as  well.  However,  a preliminaiy  numerical  study  suggested  that  the  matrix 

of  coefficients  of  the  .system  of  25  equations  has  a rank  of  17  and  consequently  only  17 

unknowns  a.  and  a.,  can  be  determined. 

> ij 

It  is  more  convenient  to  express  the  wave  drag  equation  in  terms  of  the  reduced 
variables  z..  Since  the  reduced  variables  alw'ays  take  the  levels  0,  1,  -2,  the  least 

square  fitting  procedure  once  developed  can  be  tipplied  to  any  other  values  of  the 
geometric  variables  without  modification.  .Also,  since  all  reduced  variables  are  of 
the  same  order  of  magnitude,  the  coefficient  of  each  term  m the  wave  drag  equation 
is  a measure  of  the  contribution  of  that  term  to  the  total  wave  drag.  In  terms  of  the 
reduced  variables  the  original  wave  drag  equation  takes  the  form 

(j  5 5 (1.3) 

Cr.,,,  ^ b + V b.z.  + V V b..z.z. 

DW  o Z-i  1 1 2^  ij  1 ] 

i=l  j=i  i=2 

When  the  coefficients  b^  and  b„  were  evaluated  using  the  calculated  wave  drag  co- 
efficients of  the  25  cells  and  the  baseline  in  a least-square  procedure,  Eq.  (1.3)  was 
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usi'il  tor  sor.u'  |)ri'!imi!':i r's'  "I  muiimiim  t irito^unititjns.  <»no  ol 

the  I o'ctlicU'c!  t'oiili^ii  rtilioi,,-.  u;is  otit.iiru  il  win-n  tlic  tiisol;i;^f  Ii’nuih  v.-.is  sol  to  that  ol  ttie 
baseline'  and  no  ollior  const  ra  n. Is  v.ito  imp*s«-<l.  I'ho  prialicli'il  cn-ss  section  ;it  I'.S. 

■)0  is  shown  in  l i^uri'  a In  a da-.ii-doi  lino,  Al.-.o  sho.vii  :.ic  tlie  railar  dish  by  a dashed 
line'  and  the  liaselin.e  iti'Ss  s«a  iio:.  in  a ihii"  solid  line  While  it  was  interestin'^'  U>  note 
that  the'  radonu'  requiri-ment  wa-^  !« a-uiit- iii.av  satis'u  I.  it  was  suspected  diat  tile  wron^ 
tri'iid  was  pri'dieted  since  a ri-<!ui  i!  .n  in  ilv  shape  la<.t«;r  t)  ' .ould  re*luee  the  cross- 
seetional  areas  all  the  wa\  Iroai  t!u'  nose  t<*  1 2:’.n.  Subsequent  veritieation  I'uns  by 

till-  meth(Kl  (d  c'h;i  raeteristies  program  shoiu'd  that  t!ie  wave  dra^  e<)ua[ion  was  under- 
prediclin'4  lh<'  wave  dra-^  and  tiu-  ptadieied  trend  in  b was  ineorreet.  ihis  poor  pre- 
diction was  attributc'd  to  the  rather  wide  rain^es  td  tlie  variables.  In  ortler  to  increase 
;he  prediction  aeeuraei'.  tite  three  ealeulaied  wave  drtru.  eoelTicients  by  the  verifien- 
tion  iims  and  their  corresponding  levels  ol  the  reduced  varitmles  (cells  27,  2s.  and  2d 
of  T ible  1 1 wi're  used  to  provide  thri'c  more  e(|uations  tor  tin.-  deterniinatkon  of  the 

ciH'ltieients  b.  and  b.  ..  whicli  are  tabulated  in  I' ibh'  a.  I- t-oin  Table  a it  is  si*en  that  in 
‘ 1.1 

jteneral  b is  one  order  of  maynituiie  "realer  than  b.  wldeh  is  one  onler  of  majoiitude 
greater  than  b. ..  I q.  (bit  with  the  iinprovtd  b^  iiul  to.  was  :i>tain  used  to  predict  the 
minimum  wave  dra^  I'onliaq  ration  when  the  lose  lave  length  was  set  to  that  of  the  base- 
line ami  the  radome  constramt  was  impos«'d.  The  residiio'.;  cross  section  I'f  the  mini- 
mum wave  dra"  contiau  ration  at  I'.S.  70  shown  m 1 lyiire  !•  I)v  a ht  ayv  solid  line, 

Uolh  th<’  predicted  trend  ami  \iave  d.aiv  e<K'tlk-ienl  wi'ix'  leribi-il  by  Hie  metluKl  of 
eharacteristii's  program,  as  presented  In  tlu-  r.e\t  section. 

A new  eonci'pt  has  tin  re  to  re  Ik-cii  ii't  r>  ilmed  to  tne  optimi/ati.>ti  nroeeiiure  usin>!: 
I.atin  Sf)uare  .sampling.  I Ins  l•l•n^ept  enables  the  prediction  e(iun!i(in  lo  "Icani"  from 
evjierienee  in  order  to  iinprove  its  pri'dK-.ion.  After  eai  h or  several  verifiention 
eak  m.dlofis  o|  the  predicted  eontli'iiral ions,  the  ''esultini;  ('<|uatlons  eorfespotulln}!  to 


TABLE  5.  COEFFICIENTS  OF  WAVE  DRAG  EQUATION 


b = 6.677233  X 10 
o 

VALUE  OF  b.  X 10^ 


i 

1 

2 

3 

4 

3 

6 

o 

X 

-3.96051S 

2.010586 

0.035260 

4.537457 

-0.191512 

2.672163 

VALUE  OF  b. . X 10^ 
‘J 


ESI 

2 

3 

4 

5 

2 

2.450636 

-6.113357 

7.665063 

-0.727027 

3 

0,134838 

-1.613047 

-2. 353471 

i 

j 

5.563432 

1.405971 

imiaiB 

5 


4.226010 


Best  Available  Co: 

the  v erification  rims  are  included  in  the  least-square  procedure  to  obtain  improved 
coetlicients  b.  and  b.^  of  the  prediction  equation.  In  general  one  or  two  applications 
of  this  process  should  yield  satisfactory  results.  The  improved  prediction  for  one  set 
of  constraints  may  or  may  not  lead  to  improvement  for  other  sets  of  constraints.  How- 
ever. the  same  process  cmi  be  applied  wherever  needed. 

An  additional  advantage  of  this  new  concept  is  that  nonlinear  terms  involving 
anti  /^.  can  also  lx;  taken  into  account.  With  each  addition  of  an  equation  due  to  a verifi- 
cation, the  rank  of  the  matrix  of  coefficients  usually  increases  by  one.  That  means  one 
more  b^^  of  a nonlinear  term  involving  Zj  or  can  be  included.  This  process  can  be 
repeated,  if  necessary,  until  all  quadric  terms  are  accounted  for. 
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5-  ClKOMETRiC  CONSTHAINTS  AND  WAVE  DRAG  REDUCTION 

The  wave  drag  equation  can  now  be  investigated  to  reditte  the  wave  drag.  In  the 
process  of  determining  the  minimum  wave  drag  configuration,  certain  geometric 
constraints,  such  as  minimum  cross-sectional  area  or  m:ucimum  ftiselage  length, 
must  be  satisfied.  For  each  set  of  geometric  constraints  there  exists  a minimum  wave 
drag  configuration.  In  this  section  the  optimization  procedure  and  some  minimum 
wave  drag  configurations  are  presented  and  discussed. 

a.  Optimization  Procedure  and  Geometric  Ccatstraints 

The  simplest  and  surest  way  to  find  the  minimum  wave  drag  configuration 
subject  to  a given  set  of  geometric  constraints  is  to  use  Kq.  (13)  to  calculate  the  wave 
drag  coefficients  for  all  allowable  sets  of  levels  of  the  variables  and  pick  the  set  that 

gives  the  least  wave  drag.  The  optimization  procedure  ccnisists  of  the  following 

♦ 

steps. 

(1)  Make  a numerical  search  through  the  ranges  of  all  six  variables,  using  the 
wave  drag  equation  to  calculate  the  wave  drag  coefficients  for  those  sets  of 
levels  that  satisfy  the  constraints. 

(2)  Identify  the  set  of  levels  of  ihe  variables  that  yields  the  least  wave  drag. 

(3)  Prepare  body  description  input  data  for  the  minimum  ».avc  drag  configuration. 

(4)  Use  the  Three-Dimensional  .Method  of  Characteristics  program  to  verify  the 
wave  drag  coefficient  predicted  by  tq.  (1,1). 

The  coefficients  b.  and  b.^  of  Eq.  (13)  were  obtained  by  applying  a least  square 
procedure  to  fit  the  hirpersurface  represented  by  Eq.  (13)  over  the  25  wave  drag  co- 
efficients given  in  Table  .1  and  the  first  four  coefficients  of  Table  4.  The  standard 
deviation  of  the  fit  was  about  1.02  x 10  which  corresponds  to  an  error  of  about  1.5T 


‘ The  Numerical  Search  Procedure  is  presented  In  Appendix  D. 


bast'd  on  the  mean  wave  dra>^  coefficient.  The  maximum  error  of  the  fit  over  the  29  data 
poults  was  about  2.4',(  • Hence  step  5 of  the  optimization  procedure  is  necessary  for 
wave  drag  reduction  of  4-5'  ,’  or  less  as  predicted  by  Kq.  (13). 

Four  tvpes  of  geometric  const ramts  are  considered  in  order  to  satisfy  most 
design  needs.  These  are: 

(1)  One  or  more  of  the  six  geometric  variables  take  assigned  values  within  the 
ranges  of  variation.  For  instance,  the  length  can  be  kept  the  same  as  that  of 
the  baseline  or  the  width  can  take  90^  of  the  baseline  wi.llh. 

(2)  At  one  or  more  fuselage  stations  the  configuration  cross  sections  contain 
geometric  curves  prescrilied  by  tabulated  data  of  widths  versus  elevations. 
This  type  of  constraint  is  useful  for  the  placement  of  particular  equip- 
ment, for  instance,  the  radome. 

(3)  At  one  or  more  fuselage  stations  the  configuration  cross  sections  satisfy 
minimum  area  requirements.  The  minimum  area  can  be  either  a given  area 
or  a certain  percentage  of  the  baseline  cross-sectional  area. 

(4)  The  configuration  satisfies  a minimum  volume  requirement  between  two 
given  stations.  A given  volume  or  a certain  percentage  of  the  baseline 
volume  can  be  imposed. 


b.  Minimum  Wave  Drag  Configurations 

Several  minimum  wave  drag  configurations  subject  to  various  geometric  con- 
straints are  presented  to  illustrate  the  application  of  this  program.  These  configura- 
tions were  obtained  by  using  the  optimizatiiwi  procedure,  and  the  predicted  wave  drag 
coefficients  were  verified  by  the  Three-Dimensional  Method  of  Characteristics  (3DMoC) 
program  in  some  cases. 
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Within  the  ranges  of  variation  of  the  six  variables  (when  the  only  constraint  was 
implicitly  imposed  by  the  ranges  under  consideration),  a 34.8',T  reduction  of  the  wave 
drag  coefficient  was  predicted  by  Kq.  (13)»  The  configuration  corresponds  to  the 
longest  nose,  as  might  be  e^qjccted.  Wlien  the  length  was  set  to  be  equal  to  that  of 
the  baseline,  a 15'<'  reduction  was  predicted  by  the  equation.  The  percentage  reduc- 
tion of  the  wave  drag  in  relation  to  the  increase  of  hiselage  length  is  plotted  in 
Figure  10.  It  is  seen  that  for  every  inch  the  fuselage  is  lengthened  the  wave  drag 
equation  predicts  a wave  drag  reduction  slightly  over  1%.  Since  no  constraint  other 
than  the  nose  length  was  imposed,  these  minimum  wave  drag  ccx)  figurations  have  the 
least  cross-sectional  areas  attainable  withm  the  ranges  of  variation  of  the  remaining 
five  geometric  variables.  Consequently  these  configurations  cannot  contain  the  radar 
dish  which  is  represented  by  a dashed  circle  at  F.S.  50  in  Figures  (i  and  9.  Wlten 
the  fuselage  length  was  set  to  that  of  the  baseline  and  the  radome  constraint  was 
imposed,  the  wave  drag  equation  predicted  a 14',V  reduction  of  the  wave  drag  coefficient. 
The  cross  section  of  this  minimum  wave  drag  configuration  at  F.S.  50  is  shown  by  the 
heavy  solid  line  in  Figure  9.  It  is  seen  that  the  lower  biselage  cross  sectiem  barely 
contains  tlie  radar  dish.  This  predicted  trend  that  the  wave  drag  decreases  with  the 
shape  factor  b in  the  neighborhood  of  this  minimum  wave  drag  configuration  was 
verified  by  calculations  using  the  Three-Dimensional  Method  of  Characteristics 
program.  Calculation  of  the  wave  drag  of  this  configuration  also  resulted  in  a wave 
drag  reduction  of  12.0'^,  which  agrees  fairly  well  with  the  14"  reduction  predicted 
by  the  wave  drag  equation,  considering  that  the  standard  deviation  of  the  least  square 
fit  is  1.  S'Jr . The  front  view  of  this  minimum  wave  drag  configuration  is  presented  in 
Figure  11.  which  shows  the  contours  of  the  cross  sections  from  F.S.  20  to  F.S.  230 
at  intervals  of  10.  The  front  view  of  the  baseline  is  presented  in  Figure  12  for  com- 
parison. It  is  seen  that  this  substantial  reduction  of  the  wave  drag  was  accomplished 
at  the  expense  of  reduced  fuselage  volume. 
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Two  important  general  constraints  that  strongly  affect  the  wave  drag  are  the 
fuselage  nose  length  and  the  fuselage  volume.  Therefore,  the  optiniization  procedure 
was  used  to  determine  a number  of  minimum  wave  drag  configurations  for  various 
nose  lengths  and  fUselage  volumes.  Since  the  horizontal  displacement  of  the  maximum 
breadth  line  has  an  important  effect  on  the  fuselage  volume,  nonlinear  terms  of  the 
displacement  a and  the  nose  length  i should  be  included  in  the  wave  drag  equation  for 
accurate  preaiction  of  the  influence  of  these  two  geometric  variables  on  the  wave  drag. 
Equation  (13)  was  first  used  to  obtain  two  preliminary  minimum  wave  drag  configura- 
tions which  are  tabulated  in  Table  4 as  cells  31  and  32.  The  30th  cell  in  Table  4 cor- 
responds to  the  minimum  wave  drag  configuration  subject  to  the  radome  constraint, 
which  was  discussed  in  the  preceding  paragraph.  These  32  cells  provided  32  equa- 
tions which  allowed  23  terms  (or  23  coefficients  b^,  b^.)  to  be  included  in  an  improved 
wave  drag  equation.  All  important  nonlinear  terms  involving  a and  / were  included. 

The  final  improved  wave  drag  equation,  with  six  additional  terms,  is  expanded 
into  the  form: 


Sw  * ^o  ^ ‘*1^1  *^2^2  ^ ‘'3^3  ^ '4^4  * *^5^5  ‘^6^6 


* * '■22^2^  * * '>55^5^  * 


(14) 


" ^21^221  - *^3223-2  " '’412421  " ‘>,2Z4Z2  * >^432423 


* *^522522  *’532523  '’542524  *’622622  ^ *’642g24 


The  coefficients  of  Eq.  (14)  are  tabuiated  in  Table  6.  The  optimization  procedure  using 
the  improved  wave  drag  equation  was  then  applied  to  determine  the  dependence  of  the 
wave  drag  on  the  nose  length  and  fuselage  volume.  The  predicted  percentage  reduc- 
tion of  the  wave  drag  coefficient  is  plotted  against  the  increase  of  the  fuselage  nose 
length  in  Figure  13  for  various  fuselage  volumes  (expressed  as  % of  the  baseline 
volume).  The  two  circles  mark  the  calculated  percentage  reductions  of  the  wave  drag 
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TABLE  6.  COEFFICIENTS  OF  IMPltOVED  WAVE  DRAG  EQUATION 


Additional  forms  a<lded  to  the  original  wave  drag  equation 


usinj;  the  Three-Dimensional  Method  of  Characlerisdcs  foi-  the  two  selected  configura- 
tions. Excellent  agreentent  between  the  predicted  end  calculated  wave  drag  reduction 
is  seen.  When  the  nose  length  and  fuselage  volume  were  set  to  the  baseline  values, 
very  little  wave  drag  reduction  was  obtained  within  the  variations  of  the  six  geometric 
variables  is  selected.  This  set  of  geometric  variables  and  their  ranges  were  chosen 
to  illustrate  the  wave  drag  reduction  procedure  in  such  generality  that  a wide  range  of 
configurations  similar  to  the  baseline  was  covered  and  various  illustrations  could  be 
madCi  No  specific  effort  was  made  to  tailor  the  variables  to  the  wave  drag  reduction 
of  the  baseline  configuration.  The  baseline  wave  drag  could  conceivably  be  further 
reduced  by  some  appropriately  selected  variables.  Referring  to  Figure  13  again,  some 
rule-of-thumb  percentage  reductions  of  the  wave  drag  can  be  observed.  On  the  average 
one  inch  increase  in  the  fuselage  length  reduces  the  wave  drag  by  sli^tly  over  one  per- 
cent while  one  percent  decrease  in  the  fuselage  volume  reduces  the  wave  drag  by  about 
3/4  percent. 

The  characteristics  of  the  improved  wave  drag  equation  are  shown  in  Figures  14 
through  18.  Each  of  the  figures  shows  the  variations  of  the  wave  drag  coefficient  with 
respect  to  two  of  the  six  geometric  variables,  where  the  remaining  four  are  set  equal 
to  tlie  values  of  the  baseline.  Although  it  is  difficult  to  plot  the  present  wave  drag 
function  In  a six -dimensional  space,  these  figures  at  least  enable  one  to  isolate  and 
study  the  interacting  characteristics  between  any  two  of  the  six  varialiles.  The  effects 
of  volume  are  also  shown  by  plotting  sj'mbols  on  the  curves  to  indicate  various  fuselage 
volume  limits.  Since  the  geometric  variable  b is  different  for  different  fuselage  sta- 
tions, the  parameter  b used  in  these  figures  is  actually  the  reduced  variable  which  is 
constant  for  each  given  configuration. 

Figures  14a  through  14e  show  the  variations  of  the  drag  coefficient  with  respect 
to  the  height  of  the  lower  deck  h for  families  of  f , b,  a',  a and  f,  respectively.  Within 
the  ranges  of  variables  stu<Med,  decreases  monotonically  with  h.  Figure  14a 
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iruiu  aU's  ihat  lor  a given  voUinu',  t decreases  with  llie  length  of  the  fuselage  n<ise 
f.  Along  values  of  constant  volume.  Figures  I4h,  14c  and  14fl  show  that  when  each  of 
the  respective  variables  (b,  a'  and  a>  interacts  with  h,  there  is  a definite  combina- 
tion of  variables  which  yields  a configuration  of  minimum  wave  drag.  For  a volume 
equal  to  the  baseline  value  tV'ol/  Vu!|jj  1),  for  example.  Figure  14b  shows  that  a 
minimum  occurs  near  b h 0 which  are  equal  to  the  baseline  values.  The  occurrence 
of  the  minimum  near  the  baseline  configuration  verifies  the  results  obtained  from  the 
numerical  search  where  little  reduction  of  wave  drag  from  the  baseline  w'as  obtained 
within  the  variations  of  the  selected  geometric  variables.  Lines  of  constant  a'  and  a 
in  Figures  14c  and  14d  respectively  are  nearly  parallel,  indicating  that  the  interaction 
between  a'  and  h and  the  interaction  Ixjtween  a and  h arc  weak.  This  fact  is  also  ^ 
revealed  by  examining  the  magnitude  of  the  coefficients  of  the  wave  drag  equation  tab- 
ulated in  Table  6.  The  magnitude  of  these  interaction  terms,  although  included  in  the 
formulation  of  the  wave  drag  equation,  is  one  order  of  magnitude  smaller  than  the 
average  of  the  remaining  interaction  terms.  Figures  I5a  and  1.>I>  show  the  varia- 
tions with  respect  to  the  length  of  the  nose  t for  families  of  a'  and  f.  re sixjc lively. 

Again,  the  trend  shows  that  decreases  with  f , but  mininia  for  given  constant 
volumes  are  not  definitely  shown.  Figure  I5b  shows  that  the  width  of  the  fuselage 
bottom  flat  .section  has  little  effect  on  the  drag  eoefficient  when  the  remaining  four 
variables  are  set  equal  to  the  values  of  the  baseline.  Figures  Kia  through  IGd  show 
the  variations  with  respect  to  the  horizontal  displacement  of  the  maximum  breadth  line 
a for  families  of  1,  b,  a'  and  f,  respectively.  Here,  the  figures  show  that  the  drag 
coefficient  increases  monotonically  with  a.  Similar  to  the  results showm  In  Figures 
11a  through  14e,  for  constant  volume  the  drag  coefficient  decreases  with  f , and  when 
the  volume  i.s  equal  to  that  of  the  ba.seline,  minima  occur  at  or  near  the  baseline.  In 
Figure  H5b,  for  example,  the  minimum  occurs  at  a = l.Ol  and  I)  - lAi  approximately 
for  Vol/Volj^j  - 1.  The  curves  in  Figures  l«a,  Kic,  and  Ifid  are  parallel  since  inter- 
action terms  lx?tween  each  of  the  variables  with  a are  not  included  in  the  present  wave 
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(inifi  i.‘(inati<<n.  Similar  to  the  rosults  shown  in  ^'i^ure  IS!),  varies  little  with  t. 

l-  igiires  17a  and  171)  show  the  variations  with  respect  to  the  vertical  displacement  of 
the  maximum  breadth  line  for  families  of  h and  f,  respectively.  Here,  a minimum 
configurtition  for  a given  volume  constraint  cannot  be  determined  graphically  but  a 
minimum  can  l)c  located  for  each  constant  b <ir  f.  Figure  I7a  shows  that  the  drag 
coefficient  increases  with  b,  except  when  a'  is  Ixdow  the  value  of  the  baseline  (a'  - 1) 
and  when  b is  less  than  about  *1/3.  Figure  I7b  shows  that  the  drag  coefficient  in- 
creases with  f when  a'  is  ixtlow  the  baseline  value  (a'  = 1)  and  the  trend  reverses  when 
a'  is  above  the  baseline  value.  Figures  ISa  and  18b  show  the  variation  with  respect  to 
the  shape  factor  b for  families  of  f and  f,  respectively.  Again,  Figure  18a  shows  that 
the  drag  coefficient  decreases  with  £.  Similar  to  the  characteristics  of  vs.  a' 
shown  in  Figure  17b,  the  trend  of  the  wave  drag  coefficient  reverses  from  increasing 
with  f when  b is  below  the  baseline  value  (l»  - 0)  to  decreasing  with  f above  the  baseline 
value.  Along  the  constant  baseline  volume  line.  Figure  18b  shows  a slight  reduction 
of  the  drag  coefficient  from  the  baseline  when  f 1.5  and  b = -.25. 
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FIGURE  17b.  CHARACTERISTICS  OF  THE  WAVE  DRAG 
EQUATION  - Crx,„  VS  a'  (FOR  CONSTANT  f) 


6.  CONCLUDING  REMARKS 


The  following  remarks  and  conclusions  can  be  made  for  Phase  I of  the  wave 
drag  reduction  program. 

a.  The  optimization  procedure  using  the  Latin  Square  sampling  technique  and 
the  Three-Dimensional  McIIkxI  of  Cluiracteristics  is  useful  and  practical. 

The  procedure  can  be  used  for  other  optimization  purposes  besides  wave 
drag  reduction, 

b.  The  geometric  variables  and  their  ranges  of  variation  should  lx*  selected 
carefully  since  the  quality  of  the  procedure  and  results  are  affected  by  the 
selection.  Some  guidelines  on  selecting  the  variables  and  ranges  are  given 
in  this  appendix,  but  etperience  is  valuable. 

c.  A new  concept  has  been  introduced  to  enable  the  optimization  procedure  to 
improve  during  the  course  of  determining  the  optimum  configuration.  This 
concept  is  particularly  useM  udien  the  variables  and  their  ranges  were  not 
well  chosen  either  because  of  inexperience  or  for  the  sake  of  compromise. 

d.  Within  the  ranges  of  variation  of  the  six  chosen  geometric  variables,  the 
following  rule-of-thumb  percentage  reductions  of  the  wave  drag  are  ob- 
tained for  the  F-4  ftiselage.  For  every  inch  the  nose  is  lengthened,  the 
wave  drag  is  reduced  by  sli^tly  over  one  percent.  For  each  percent  the 
fuselage  volume  is  decreased,  the  wave  drag  is  reduced  by  about  3/4  percent 
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APPENDIX  B 


THE  LATIN  SQUARE  METHOD 

1-3 

The  Latin  Square  is  a samplinfr  method  for  designing  an  experiment. 

An  N—  order  Latin  Square  is  an  arrangement  of  the  N values  of  N 1,  or  fewer, 
variables  (where  N is  a prime  number)  in  such  a way  that  the  resulting  surface 
fitting  process  is  efficiently  performed.  The  method  may  be  used  for  determining 
an  approximate  equation  for  a function  of  a number  of  variables  and  hence  for  solv- 
ing extremum  problems  when  the  (tependent  variable  is  not  known  as  an  explicit 
function  of  the  independent  variables. 

The  dependent  variable,  in  this  case  wave  drag  coefficient  is  assumed  to 

be  a quadric  function  of  the  independent  variables  x.  (which  are  related  to  the  geometry 
of  the  forward  fuselage  in  this  application), 

6 5 5 

* 2 “1*1  * S 2 “li  *l*i 

i = 1 j = I I = 2 

Equation  (1)  is  written  for  a 5 x 5 Latin  Square  with  six  independent  variables,  i.e. 

N »•  5.; 


For  a given  set  of  independent  variables,  Xj,  a new  set  of  reduced  variables, 
say  Zj,  may  be  determined  such  that  the  variables  Zj  take  the  values  0,  ±1,  t2, 

N-1 

. . . Jtct-'Ofding  to  the  transformation  equation 


The  5 X 3 Latin  Square  arrangement  in  terms  of  the  levels  of  the  reduced 

variables  Is  shown  in  Table  i . The  arrangement  may  appear  to  be  at  once  regular 

i 


and  random.  It  appears  to  !>e  regular  since  all  variables  t;^e  orderly  permutated 
values.  It  appears  to  be  r;mdom  since  there  seems  to  be  no  bias  in  re{)resenting 
all  segments  of  the  population.  Regularity  and  randomness  - therein  lies  the  ef- 
ficiency of  the  Latin  Square  sampling  technique.  It  has  been  shown  by  Yates^ 
that  randomized  blocks  are  more  efficient  than  complete  randomization  and  Latin 
Squares  are  more  efficient  than  randomized  blocks.  He  also  found  that  on  one  set 
of  experiments  about  2 1/2  times  as  many  plots  were  needed  with  the  randomized 
blocks  of  5 blocks  as  with  the  5x5  Latin  Squares  to  obtain  the  same  accuracy  (see 
Reference  1,  pp.  202-203).  Careful  examination  of  Table  1 revealed  that  in  all 
the  cells  having  a given  level  of  any  reduced  variable  any  one  of  the  remaining 
variables  takes  on  all  possible  levels  once  and  only  once.  For  example,  if  we 
select  Zj  = -2  to  be  the  given  level  of  the  one  reduced  variable,  then  column 
one  contains  all  cells  having  Zj^  = -2.  It  is  seen  that  in  column  one  each  of  the  five 
levels  of  any  one  of  the  remaining  variables,  say  z„  appears  once  and  only  once. 
Consequently,  the  number  of  all  cells  having  a given  level  of  any  reduced  variable 
is  five  since  there  are  five  levels  for  each  of  the  other  variables  to  take  through  a 
permutation  process.  Several  consequences  can  be  observed.  Firstly,  if  we  com- 
pare any  two  cells  we  find  that  only  one  of  the  variables  has  the  same  level  in  both 
cells,  whereas  all  other  variables  have  different  levels.  Secondly,  for  each  re- 
duced variable  there  arc  five  groups  of  five  cells  which  satisfy  the  condition  that  in 
each  group  this  reduced  variable  takes  one  of  the  five  levels.  Thirdly,  by  a suitable 
rearrangement  of  the  cells  any  pair  of  the  reduced  variables  can  take  the  column  and 
row  indices  as  their  levels  just  as  Zj  and  do  as  shown  in  Table  1.  To  do  this, 
arrange  each  group  with  the  same  level  of  one  variable  as  a column  and  arrange  the 
columns  in  an  ascending  order  with  respect  to  the  levels  of  this  variable.  Within 
each  column  (or  group)  arrange  the  cells  in  an  ascending  order  with  respect  to  the 
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levels  of  the  other  variable  of  the  pair.  These  properties  slunv  heuristieally  the. 

effic  ieney  of  the  Latin  S<iua re  sainplinii  teehnique  which  has  been  used  tor  develop-  . 

4-5 

mt;  various  optimal  designs. 

It  may  seem  intricate  to  construct  a Latin  Square  of  this  type,  especially  when 
the  size  exceeds  5 x o.  The  followina  steps,  however,  '-an  be  followed  to  construct  . 
this  type  of  Latin  Square  of  any  size  as  lonp,  as  the  iiumlier  of  rows  (or  colunins)  is  a 
prime  numlx'r.  The  steps  are  best  explained  usinn  Table  1 as  an  example. 

(a)  The  variables  Zj  and  are  assigned  levels  equal  to  the  column  and  row 
indices  as  shown  in  Table  1. 

(bl  For  each  of  the  diagonal  cells  where  /.^  - z^.,  the  level  of  z.  is  increased 
successively  from  i 1 to  i (*  by  a level  equal  to  that  of  Zj  in  a cyclic 
order,  in  which  the  lowest  level  (i.c..  -2)  is  considered  one  level  higher 
than  the  highest  level  (i.e.,  2). 

(c)  The  remaining  levels  of  z.,  to  in  each  tow  are  obtained  by  permutation 
in  an  increasing  cyclic  order  using  the  diagonal  cell. 
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APPENDIX  C 


BODY  DESCRIPTION  METHOD 


In  this  Appendix,  the  general  body  description  method  is  described  in  some 
detail  and  is  illustrated  by  the  body  description  of  phase  I.  The  technique  of  fitting 
a conic -section  curve  is  also  given. 

The  fuselage  is  described  by  a number  of  generating  lines,  such  as  the  upper 
profile,  the  lower  profile,  the  maximum  half-breadth  line  and  generating  lines  for 
the  canopy.  In  the  present  body  description  procedure,  each  generating  line  is 
divided  into  as  many  segments  as  necessary  to  allow  each  segment  to  be  described 
by  a conic-section  curve.  At  each  cross  section  of  the  fuselage  simple  analytic 
curves,  e.g. , the  ellipse  or  cubic,  connect  any  two  adjacent  generating  lines  to 
form  the  contour  of  the  cross  section.  The  fuselage  is  thus  described  analytically 
by  simple  low  order  curves. 

The  body  description  procedure  detailed  here  was  developed  for  use  in  con- 
junction with  the  three-dimensional  method  of  characteristics  for  calculating  flow 
fields  over  smooth  bodies.  For  a smooth  body  a unique  normal  to  the  surface  must 
exist  every  where,  and  this  condition  of  a unique  normal  usually  requires  that  the 
curves  and  their  slopes  be  continuous  at  the  juncture  between  two  contour  curves 
or  two  segments  of  a generating  line. 

1.  DEFINITION  OF  GENERATING  LINES 

The  fuselage  is  located  in  a right-handed  coordinate  system  where  the  Y-axis 
is  aligned  with  the  fuselage  axis,  the  X-axis  is  spanwise  and  the  Z-axis  is  up.  A 
schematic  of  a fuselage-canopy  configuration  is  shown  in  Figure  1 and  a typical  cross 
section  in  Figure  2. 

In  general,  five  generating  lines,  each  being  described  by  two  functions  of  Y 
(representing  their  vertical  and  horizontal  projections  and  -Xj,  respectively,  by  a 
general  curve  fit  of  conic  section^  plus  an  additional  three  lines,  each  being 
described  by  a single  function  of  Y,  serve  to  define  the  fuselage-canopy  configuration 
in  phase  I. 


Cl 


Those  are  identified  in  Fi^iires  1 :uid  2 as 

(u)  I’pper  i)rofile  line 

(0)  Maximum  half-hreuilth  liiic 
® Lower  profile  line 

© Canopy  crown  line 

Cimopy  definition  line 
Shajio  factor  b,  loctitin^  the  45®  tangent  line 
© Upper  fairing  line 

© Lower  fairing  line 

where  lines  ®,  ®,  ©,  © arc  required  only,  of  course,  in  the  region  of  the  canopy. 
The  conic-section  curves  describing  the  lines  take  the  form 

(|)  - PY  ♦ a t IRY**  SY*T  (1) 

Details  of  the  determination  of  the  coefficients  P,  Q,  R,  S,  T of  the  general  conic 
curves  will  be  discussed  in  Section  4.  Each  conic-section  segment  must  be  con- 
tinuous with  the  previous  segment  and,  with  very  few  exceptions,  the  slope  must  be  , , 
continuous  at  the  junctures  to  satisfy  the  requirement  of  a unique  normal  to  the  sur- 
face discussed  above,  l.'nlike  others,  line  © is  not  a true  generating  line.  It  is 
actually  a shape  factor  for  the  lower  part  of  the  fuselage,  which  expresses  the  dis- 
tance to  a tangent  line  as  a function  of  Y.  A full  discussion  of  the  details  is  pre- 
sented in  Section  2(c). 

Lines  © and  © which  fix,  at  each  cross  section,  the  extent  of  the  fairing  cubic 
between  canopy  .?nd  upper  fuselage  are  defined  by  their  projections  in  the  YZ-planc, 
i.c.  Zf,(\')  and  Zfj(Y),  only.  Their  projections  and  Xf,(Y)  are  constrained 

to  lie  on  the  ellipses  defining  the  can<j|)y  and  upper  half  of  the  fuselage,  Respectively, 
as  detail c.l  in  Section  2(d). 


2.  CONTOUR  CURVES  FORMING  CROSS  SECTIONS 

As  shown  in  Figure  2,  for  the  most  general  case  of  a fuselage-canopy  configunition 
with  a windshield  flat,  the  contour  of  the  cross  section  begins  with  a striUght  lino 
representing  the  flat  from  the  centerline  to  ©.  The  canopy  contour  from  ©to  © 
is  circular,  but  can  be  elliptic  in  general.  The  upiier  fuselage  is  represented  by  an 
elliptic  curve  between  (m)  and  ©.  The  lower  fuselage  is  represented  by;a  general 
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conic-scction  curve  from  0 to  0;  the  bottom  flat  by  a sti-aight  line  from  0 to  the 
ccmerlinc.  The  intersection  between  the  cano|>y  and  the  fuselage  is  faired  by  a cubic 
from  0 to  ©. 

a.  Canom’ 

The  cross-sectional  shape  of  the  cant^y  is  assumed  to  be  circular.  If  there  is  no 
windshield  flat,  0 will  lie  on  the  centerline  (Figure  3a)  and  the  equation  for  the 
cmiopy  curve  is  given  by 

‘ -0 

If  a windshield  flat  does  exist,  0 moves  off  the  cafiterline  as  shown  in  Figure  3b  and 
Equation  (2),  describing  the  canopy  between  0 and  0 leads  to  a slope  discontinuity 
at  0 , the  juncture  between  the  circular  arc  and  the  windshield  flat.  Thus,  a quarter 
ellipse  is  used  to  approximate  the  circular  arc  between  0 and  0 to  give  a slope  of 
dZ/dX  - 0 at  0,  providing  slope  continuity  with  the  flat  portion  of  the  contour.  The 
equation  for  this  approximating  ellipse  is 


The  details  of  determining  the  point  0,  given  the  locations  of  the  canopy  crown  and 
sill  lines,  are  presented  in  Section  3. 

b.  Upper  Fuselage 

The  equation  for  the  upper  fuselage  is 

indicating  that,  for  the  Ix>dy  description  in  its  present  form,  X is  constrained  to  be 
zero  for  all  Y. 

Again  it  may  be  noted  that  0,  in  the  region  of  the  canopy,  is  fictitious  and  cannot  be 
obtained  directly  from  three-view  drawings.  Details  of  determining  0,  given  the 
location  of  the  maximum  half-breadth  line  0 and  the  location  of  the  canopy  sill  line, 
may  be  found  in  Section  3. 
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c.  Lower  Fuselage 

Referring  to  Figure  2,  the  contour  curve  from  @ to  has  a zero  slope  at  © and 
the  slope  approaches  infinity  at  As  in  aircraft  lofting  practice,  the  shape  of  this 
contour  curve  is  determined  by  specifying  the  distance  from  the  origin  to  one  of  its 
tangent  lines  that  makes  a 45“  angle  with  the  X-axis.  Hence,  this  distance,  which  is 
designated  by  b in  the  following  derivation  of  the  equation  for  the  contour  curve,  may 
be  regarded  as  a shape  factor  for  this  curve. 

The  equation  for  the  45“  tangent  line  is 

Z • X - y?  b 

which,  wh«j  the  origin  is  translated  to  0t  becomes 

Z'  - X'4-  H (5) 

where  Z'  = Z-Z„„  X'  = X-X^  and  H * Xm-  Zm-  ^/2  b.  In  the  new  coordinate  system,  the 
the  general  quadric  equation,  satisfying  the  conditions  X'(Z^)  * X^,  Z'(0)  = 0, 

(dX'/dZ')m  =0  and  (aZ7aX')^  »0,  reduces  to 

K(X/Z'- Z)X'>*  ♦ X'{Z'-ip  - 0 (6) 

which  represents  a family  of  curves  with  K as  a parameter.  To  determine  K,  Eq.  (5) 
is  substituted  into  Eq.  (6)  to  5deld  a quadratic  equation  of  the  form 

Ax'*+8x'+C  - 0 (7) 

The  condition  for  Equation  (6)  to  be  a tangent  line  is  that  Equation  (7)  have  a double 
root:  i.  e.  B’-4AC  » 0,  which  leads  to 

K . 

Equation  (6)  with  K given  by  Eq.  (8),  then  represents  the  contour  curve  from  @ to  0. 
The  range  of  variation  of  this  curve  obtainable  by  applying  this  equation  is  illustrated 
in  Figure  4,  where  a family  of  conic  section  curves  is  given  for  different  values  of  b. 

By  increasing  the  distance  b,  the  curve  is  seen  to  vary  from  almost  a straight  line 
to  a sharply  bending  curve  approaching  the  two  sides  of  a right  triangle. 


4^X/Z/(M*X/-Zy. 
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d.  Fairing  Curve 


Thi-  intersection  between  the  canopy  and  upper  fuselage  is  faired  by  a cubic  from  (|i) 
to  ©.  The  projections  Zf,(Y)  and  Zf,(Y)  of  lines  © and  © on  the  YZ-plane  are  given 
by  conic-section  curves  (such  as  Equation  (1),  where  P,  Q,  R,  S,  T are  input  quantities). 
The  projections  Xf,(Y)  and  Xf,(Y)  on  the  XY-plane  are  obtained  by  solving  Eqs.  (3)  and 
(4),  respectively. 


(9) 

The  fairing  curve  matches  the  slopes  of  the  ellipses  at  both  end  points  fl  and  f2.  The 
slopes  are  defined  by  differentiation  of  Eqs.  (3)  and  (4) 


a*  VXri-XcMXc-Z*^ 

2£jei<y)  . - (IcLllaW 

02  ' Xft  ' ' 


(10) 


The  cubic  equation  that  satisfies  the  conditions  of  Eqs.  (9)  and  (10)  at  point  ft  can  be 
written 


X-  X*.»  OXfl(Z.Z>)  ♦ d(2“2fi\*»el  (Z’Zfi)*  (11) 

The  coefficients  c and  d are  obtained  by  applying  the  conditions  of  Eqs.  (9)  and  (10) 
at  point  f2 


C-  3(Xu-Xh)  *2  ^)(Zu-Zfl)  ) 

(12) 

d - -2(X„-Xu)+ (2|l*  ^)(Z,t-2q) 

Equation  (11)  with  c and  d given  by  Eq.  (12)  is  then  the  cubic  equation  for  the  fairing 
curve  from  © to  ©.  The  quantities  Xf,,  Xf,,  0Xf,/0Z,  OXfj/OZ  are  given  by  Eqs.  (9) 
and  (10),  while  X^  Xg,  X^,  Z„,  Z,^,  Z^,  Zf,,  and  Zf,  are  obtained  from  the  genera- 
ting lines  input  in  the  form  of  Equation  (1). 
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3.  DETERMINING  GENERATING  LINES  AS  FUNCTIONS  OF  Y 

The  generating  lines  defined  in  Section  1 are  determined  primarily  by  utilizing  top 
and  side  view  drawings  as  well  as  available  cross-sectional  views.  The  top  and 
side  views  arc  used  to  obtain  the  variation  of  generating  lines  (u)  and  (|)  in  the 
X Z-plane  and  of  line  @ in  the  XY-p!ane.  The  variation  of  line  (m)  in  the  YZ-plane 
and  of  line  Q)  in  the  XY-planc  are  obtained  primarily  by  utilizing  cross-sectional 
views,  as  is  the  description  of  the  shai>c  iiarametcr  ® as  a function  of  Y. 

From  Figure  2 it  can  be  seen  that  line  ®,  in  the  region  of  the  canopy,  and  line  ® 
are,  in  effect  fictitious  lines  whose  description  three-view  drawings  will  not  yield  di- 
rectly. The  two  lines  are  determined  by  assuming  that  the  cross  section  of  the  canopy 
is  a circular  arc  and  by  constraining  the  intersection  of  this  circular  arc  and  the  ellipse 
defining  the  upper  portion  of  the  fuselage  to  lie  on  the  canopy  sill  line,  which  in  turn  is 
well  defined  from  available  three-view  drawings. 

In  Figure  3,  let  Xj,  2^  be  the  location  of  the  canopy  sill  line  for  the  cross  section 
shown.  Utilizing  the  geneial  equation  for  an  ellipse  with  center  at  Xq*  Zq  and  semi- 
major and  minor  diameters  a and  b,  respectively. 


one  may  determine  Zy  by  finding  b in  the  above  equation  for  the  ellipse  through  the 
point  Xj,  Zj  with  known  center  X©  *=  0,  Z©  ^ Zn,  and  a = Xn»,  yielding 

Zu  - Zm*  (13) 

To  determine  Z^  and  X^,  defining  one  end  point  of  the  quarter  circle  passing  through 
Zc,  Xg  and  Zj,  Xi  whose  center  lies  on  Xq  - 0.  it  is  necessary  to  obtain  Zq  and  r 
for  the  circle  satisfying  the  following  conditions 

X*  ♦ (Zi-Z.)*  <•  r*  (14a) 

Xj  ♦ (Z«-Z.)*  - r*  (14b) 

Solving  for  Zg  yields 

Z^  . 2o  , Xc^^zSZv  (15) 

2(Zc-Zi') 


CIO 


and  using  Equation  (14b)  to  solve  for  r,  gives 

X,  - r - I Xc*  (16) 

Whi'n  a windshield  or  canopy  flat  exists  (i.  e.  Xc  *0,  as  shown  in  Figui*e  3b)  a slope 
disc  ontinuity  results  at  X^,  Z^.  Thus,  for  ^0,  the  circular  arc  is  replaced  by 
a quarter-ellipse,  passing  through  X|,  Zj  and  one  end  point  being  Xg,  Zc  with  dX/dZ  =0 
to  I'rovide  slope  continuity. 

This  ellipse  is  intended  to  be  an  approximation  to  the  circular  arc  determined  previous- 
ously.  It  is  obtained  by  decreasing  the  canopy  flat  (l.e.  X^)  by  an  arbitrary  percen- 
tage and  solving  for  a new  Z^  and  X^  by  an  iterative  procedure,  utilizing  Z^  and  X^ 
of  Eqs.  (15)  and  (16)  as  a first  approximation,  while  constraining  the  product  of  the 
two  semi-axes  of  the  ellipse  to  be  equal  to  r*  of  the  initial  circular  arc. 

Given  a tabulation  of  X^,  Z^,  and  X.,  Z.,  values  of  Z^,  Z ^ and  X ^ can  be 

determined  for  a number  of  faselage  stations;  their  fit  by  a sequence  of  conic-section 
curves  as  a function  of  Y can  then  be  accomplished  in  a straightforward  manner  by 
the  techniques  described  in  the  following  section. 

4.  EVALUATION  OF  THE  COEFFICIENTS  FOR  A CONIC  SECTION 

The  general  equation  of  the  second  degree,  which  represents  a conic  section,  is 

ay*  ♦ bxy  ♦ cx*  ♦okj+cx+J’*0 

utilizing  the  coordinate  system  of  Figure  5.  The  general  equation  appears  to  have  six 
constants,  but  in  reality  has  only  five  arbitrary  constants,  since  Equation  (17)  can  be 
divided  by  any  of  the  constants,  leaving  an  equation  with  five  undetermined  coefficients. 
Therefore  five  independent  conditions  are  sufficient  to  determine  a conic. 

In  the  following  sections,  two  ways  In  which  these  five  independent  conditions  are 
supplied  will  be  treated.  Section  4(a)  deals  with  the  conic  determined  by  five  }X)ints, 
no  three  of  which  are  collincar.  Section  4(b)  treats  the  conic  determined  by  one  point 
and  two  point  slopes.  By  a |x>int  slope  it  is  meant  that  a point  is  given  and  the  slope 
of  the  curve  at  the  point  is  also  given. 


a.  Conic  Determined  by  Five  Points 


Consider  the  ])roblem  ol  lindin*;  the  etjuation  of  a conic  through  five  given  points,  no 
three  of  which  are  collinoar.  Let  the  joints  he  1,  2,  :>,  4,  5 as  shown  in  Figure  5a. 


Then, 

y - miX  - hi  » 0 j 

y - rrtaX  - hi  - 0 ^ {IS) 

y - m^x  - hj  - 0 ( 

y - m4  X - h^  • 0 

represent  the  equations  of  the  lines  connecting  points  1-4,  1-2,  2-3,  3-4,  respectively. 
Then  the  equation 


Xfy-m,x -h,Ky -mjX -hj)  + (y-mjX-hjKy-m^x - o (19) 

represents  the  family  of  conics  through  the  points  I,  2,  3,  4.  In  this  equation,  K is  a 
parameter  whose  value  can  be  determined  from  the  fifth  condition,  which  is  that  the 
conic  pass  through  point  5.  (If  point  1 is  substituted  into  Equation  (19),  the  result  is 
identically  equal  to  zero  since  point  1 was  used  to  evaluate  the  coefficients  m,,  h,  and 
mj,  hj.  The  same  argument  holds  for  points  2,  3,  4.) 

K « -hiXys -m^Xs  “h4>  (20) 

” (y5“"»jX5  ■*h,Ky5-m5X5 -h,) 

After  expanding  Equation  (19)  and  collecting  like  powers,  comparison  with  Equation  (17) 
shows  that 

a m K * t 

b«  • *m4  ] 

c « 

d « - lK(ht*  hj)  w hi  ♦ h<  1 
a • K(mthi*  mshi)  + 
f m Khihi*hih4 

Rewriting  Equation  (17)  in  the  form 

ay*  (fax  ♦ d)  y ♦ cx*  + ex  + f ••  0 

and  solving,  using  the  quadratic  formula,  yields 

y . - (tax»d)  i [(bx^d^.  4a(cx 

Za 


C12 


' 1 


y 


_ Jj.x  -JSL  ♦ f ♦ 2 bdx  *■  _ (4acx*-»  4aax  * 4af~) 

2a  2a  “ I 4a*  4i^* 

b X f b^x^-*- 2bdx -Ki^  _ CX*  » CX  4 f 1 
2a  2a  " I 4a*  a * 


-Jz.X 

2a 


/ ^ _ £.\x 
V 2a*  a) 


(22) 


whifh  can  be  recognized  as  Equation  (1),  with 

pm-  b/2a 
Q m - d/2a 
R - P*-  c/a 
S - - Pd/a  - «/a 
T - Q*-  f/a 

and  the  sign  of  the  radical  to  be  determined  by  evaluation. 


b.  Conic  Determined  by  a Point  and  Two  Point  Slopes 

Now  consider  the  problem  of  finding  the  coefficients  of  a conic  determined  by  one  point 
and  two  point  slopes.  This  establishes  a curve  that  is  tangent  to  a given  line  at  a given 
point,  is  tangent  to  another  line  at  a given  point  and  passes  through  a third  given  point. 
This  particular  manner  of  determining  a conic  is  especially  suitable  since  it  allows 
for  specification  of  slope  continuityat  the  two  end  points  of  the  conic-section  segment, 
wnile  the  great  latitude  afforded  by  the  third  point  (control  point)  in  fixing  the  shape 
of  the  curve  is  utilized  to  good  advantage  in  obtaining  a good  aj5i>roximation  to  the  data 
to  be  fitted. 

Compare  Figures  5a  and  5b.  The  )K>ints  1 and  2 of  Figure  5a  coincide  in  Figure  5b, 
and  the  points  3 and  4 of  Figure  5a  coincide  in  Figure  5b.  The  line  which  joins  jxjints 
2,  3 and  the  line  which  joins  points  4,  I in  Figure  5a  therefore  coincide  in  Figure  5b. 
The  lines  1-2,  and  3-4  of  Figure  5a  become  tangents  to  the  conic  in  Figure  5b.  Thus, 
it  can  i)e  seen  that  the  .special  conditions  a|»plying  to  Equation  (21)  are  m,  = mj  and 
h,  hj,  and  the  constants  for  the  conic  may  be  written  as 

a m K * i 

b - -UKrtti  * nil  * nt4  ] 
c • Knti  * 

d - *l2Kht  * hi  * ft4  J 

a m * ntjh,  * m^ht 

f m Kh\  * hzh4 
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whi  re  the  value  of  K is  now  given  by 


X « . fyg-fHzXs  -HaKy,-  m^Xs  -h^)  (25) 

fys-m<xV-^i)* 

The  coefficients  P,  Q,  R,  S,  T for  the  conic  in  form  of  Equation  (1)  may  now  be 
coirputed  by  using  Eq.  (25)  and  Ecj.  (24)  in  Eq.  (23). 

c.  Special  Forms 

Two  special  forms  of  Equation  (1)  are  worth  noting  due  to  their  frequent  utilization 
in  this  body  description  procedure. 

A straight  line  has  the  coefficients 

P a alopc 
Q - cortftUinl 
R • 3 » T • 0 

A circle  of  radius  r«,  with  its  center  at  x,,  y*,  has  as  its  coefficients 

P « 0 
Q - yo 
R - -1 
S - 2Xo 

T - 
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APPENDIX  D 


NUMERICAL  SEARCH  FOR  MINIMUM  PROCEDURE 


The  numerical  search  procedure  to  determine  the  minimum  wave  draj;  configura- 
tion involves  the  computation  of  wave  drag  co<  fficients  of  a large  number  of  discrete 
points  blanketing  the  entire  region  in  the  six-dimensional  space  bounded  by  the  given 
ranges  of  variables.  The  point  which  has  the  lease  wave  drag  is  taken  to  be  the  mini- 
mum wave  drag  configuration.  The  accuracy  of  locating  the  minimum  is  th>  -eforc 
dependent  on  the  resoltuion  of  the  numerical  network.  Accurate  location  of  the  mini- 
mum can  be  obtained  efficiently  by  a search-by-steps  technique  where  the  search  is 
first  conducted  using  coarser  grids  to  determine  an  approximate  location,  and  then 
repeated  in  a smaller  region  centered  around  the  location  of  the  previously  obtained 
minimum.  This  process  ca.i  be  repeated  as  many  times  as  needed  to  obtain  the 
desired  accuracy.  The  major  advantage  of  such  a numerical  search  procedure  is  that 
the  given  constraints  can  be  satisfied  easily  tn’  merely  rejecting  points  which  violate 
the  constraints  during  the  search. 

The  searching  procedure  is  summarized  as  follows: 

1.  Select  the  ranges  of  search  by  specifying  the  lower  and  upper  limits  of 


each  of  the  six  reduced  variables; 


j i.« 


Select  the  resolution  for  each  of  the  vari:ibles  by  specifying  the  numbers 


of  node  points  Nb  (j  1,«)  of  that  variable.  The  rcsol'ii")'i  is  given  Ijy 


resolution 


Nj.l 
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l lii*  mimlu  r of  points  for  cnc-li  sonrrli  (.vi  lr  is  rqiuil  to  tlic  prixluct  of 


N.,  i.i' 


U 

total  points  = Xj 


For  oxampli',  if  10  points  arr  siK-cifird  fr>r  all  variables,  the  total  numixtr  of  points 
is  lo'\ 

;t,  t'ottipute  the  six  taldes  of  variables: 


i ..Nj 


('ll.  ('LljM's),  (/•„), 
4.  At  every  possible  eombination  of  the  six  vari;il)les, 

C'O,  •(''), 


j * 1,  6 


evaluate  the  given  constraints  (e.  g. , the  volume).  If  the  combination  satisfies  the 
constraint,  ct)mpute  the  wave  drag  coefficient  using  the  wave  drag  equation. 


5.  Determine  the  set  of  variables  which  yields  the  least 
efficient.  The  configuration  defined  by  z'j,  z'^,  . , . 


wave  drag  co- 


z'  is  the  mini- 
() 


mum  wave  drag  configuration  unless  one  wishes  to  improve  the  accuracy  by 
repeating  the  search.  ! 

I 

6.  If  the  search  is  to  lx*  repeated,  select  the  new  ranges  of  the  six  variables 
by  redefining  the  lower  and  upper  limits  as, 


’N.  ')  ■‘j 


where  d.  are  arbitrary  increments  of  ZjSelected  to  define  a smaller 

* J i 

region  of  search.  In  situations  where  f \ or  l z \ defined  by  the  above 


relations  are  lx.'vond  the  ranges  of  consideration,  they  are  then  set  equal 
to  the  minimum  r maximum  aceordinKly. 

Itepoat  the  prfx?cdure  from  step  (2)  until  the  desired  accuracy  is  obtained. 


APPENDIX  E 


DEMONSTR.\TION  OF  LATIN  SQUARE  PRIXEDURE 


An  IR&D  program  has  been  carried  out  to  demonstrate  the  Latin  Square 
optimization  procedure  using  the  von  Karman  ogive  as  the  baseline  configuration. 
The  area  distribution  of  the  slender  bodies,  which  yield  the  ogive,  can  be  expressed 
as 


where  i is  the  length,  B the  base  area  and  9 = » represents  the  nose,  9 = 0 the 

base;  the  a_’s  are  the  Fourier  coefficients  of  the  lineal  source  distribution  strength, 
n 

The  wave  drag  is  given  by 


i.  p 
2 


JL 

ui 


+ 2a„  + 3a«  + 4a^ 
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The  Latin  Square  procedure  is-  demonstrated  in  the  following  way:  Select  a set 
of  geometric  variables  and  their  ranges  of  variation.  This  selection  defines  a 
family  of  configurations,  over  which  the  optimization  procedure  is  then  applied 
to  yield  the  configuration  with  the  least  wave  drag.  If  the  procedure  works,  the 
minimum  wave  drag  configuratimi  will  be  a good  approximation  to  the  von  Karman 
ogive  obtained  analytically. 


El 


For  the  5x5  Latin  Square,  the  six  variables  and  their  ranges  are: 


.962  s 1 « 1.038 

.055  tf  B ^ .079 


.02  < S .02 


.01  ■?  s .01 

.02  s a.  $ .02 

4 

.02  ^ a^  S .02 


The  family  of  25  configurations  defined  by  these  variables  is  shown  in  Figures  1 
and  2,  where  the  dotted  lines  represent  the  base  line  at  I = 1.0  and  B = 0.067, 
which  were  selected  to  reflect  the  variables  used  for  the  F-4  during  phase  I studj'. 

The  demonstration  of  the  Latin  Square  procedure  and  our  finds  are  sum- 
marized below.  A straightforward  application  of  the  original  Latin  Square  pro- 
cedure did  net  yield  the  correct  minimum  wave  drag  body.  Under  the  constraint 
of  given  length  1.0)  and  base  (B  = 0.067),  the  original  procedure  predicted 
the  minimum  wave  drag  body  to  correspond  to  ag  * 02,  a^  » . 02,  002 

and  Sg  = -.  02  while  for  the  von  Karman  ogive  all  the  a*s  vanish  *=  a^  = = 0 . 

Of  these  four  variables  only  a^  is  barely  acceptable,  considering  the  numerical 
accuracy.  However,  when  the  iterative  type  Improvemrait,  which  consists  of 
adding  more  nonlinear  terms  in  the  wave  drag  function  (see  Reference  1 of  the 
main  text),  was  applied,  the  optimization  procedure  turned  out  successfully. 

Indeed,  when  only  one  case  of  ag  ® 0,  a^  = o.  2,  a^  * 0 and  ag  ® 0 was  added  to 
yield  a square  term  for  a^  in  the  wave  drag  function,  the  improved  procedure 
predicted  the  minimum  wave  drag  body  to  correspond  to  a^  ~ -.00026,  -.  00021, 

a^  = 00006,  ag  * .00031,  which  is  the  von  Karman  ogive  recovered  numerically. 


Similar  results  were  obtained  for  diffeirent  constraints  of  1 and  B. 


'I 


In  sum,' the  improved  Latin  Square  optimization  procedure  has  numerically 


predicted  the  minimum  wave  drag  body  that  was  obtained  analytically. 
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